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Abstract  of  Dissertation  Presented  to  tlie  Graduate  Council 
in  Partial  Fulfillment  of  the  Requirements  for  the  Degree  of 

Doctor  of  Philosophy 

ENSEMBLE  Qt^R'^CrERISTICS  OF  \1SU.U 
EVOKED  CORTICU  POTENTIALS  IN  NOISE 

By 

Timothy  C.  Doyle 

August,  1969 

Chairman:  Dr.  D.G.  Childers 

Major  Department:  Electrical  Engineering 

To  date,  most  investigations  of  evoked  potentials  monitored  at 
the  human  scalp  have  utilized  the  average  of  a  series  of  responses 
recorded  from  a  single  subject  to  define  components  of  the  evoked 
response,  ER.  This  technique  has  provided  a  fairly  consistent  estimate 
of  an  individual's  average  ER.  However,  with  this  method  it  is  not 
uncommon  to  obtain  quite  dissimilar  averages  from  different  subjects 
under  identical  experimental  conditions.  This  contrast  in  the  ER 
estimates  across  a  population  has  hindered  the  interpretation  of  the 
effects  of  stimulus  parameters  on  the  ER  and  frustrated  the  attenpt 
to  use  this  average  for  clinical  diamosis. 

The  present  stud>'  investigated  one  metliod  of  obtaining  infonna- 
tion  about  the  diaracter  of  the  response  to  visual  stimuli  common  to 
a  considered  normal  subject  population.  The  estimation  technique 
utilized  was  to  average  across  an  ensemble  of  fift>'  subject  EEC's 
containing  evoked  activity.  This  method  permitted  estimation  of 
common  inter-subject  response  activity  evoked  by  a  single  stimulus. 
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IX 


By  averaging  tiie  resultant  waveforms  of  tiiis  enseinble  average  across 
stimuli,   a  reliable,   relatively  noise-free  inter-subject  response 
template  was  obtained. 

The  resemblance  of  the  ensemble  of  cortical  potentials  to  a 
classical  stochastic  process  enabled  the  author  to  investigate 
three  prevalent  assumptions  about  the  statistical  character  of  the 
evoked  response  and  background  noise  activity.     These  three  assumptions 
are:    (1)   the  response  and  background  noise  are  independent;    (2)   the 
background  noise  is  stationary ;   and  (3)   the  response  and  noise  are 
additive.     Tlie  results  presented  in  this  dissertation  lend  some 
measure  of  support  to  each  of  the  three  assumptions. 


CHAPTCR  1 


INTRODUCTION 


The  time-varying  electroencephalographic  potential  specifically 
elicited  by  a  single  controlled  sensory  stimulus  is  referred  to  as 
an  evoked  potential  or  evoked  response  (ER) .  An  invariant  evoked 
response  to  controllable  stimulus  parameters,  monitored  from  the 
human  scalp,  would  provide  an  excellent  empirical  diagnostic  tool 
for  brain  researdi.  However,  this  has  eluded  all  investigators 
because  of  two  inherent  limitations.  First,  the  parameters  affecting 
the  ER,  such  as  physiological  variability  of  the  subject  during 
stimulation,  are  usually  unkno\vn  and  uncontrollable.  Second,  the 
obsen^ability  of  a  single  response  is  severely  handicapped  by  the 
interaction  of  the  desired  signal  (ER)  with  the  background  bioelectric 
on-going  activity. 

The  physiological  variability,  or  state  of  the  subject,  intro- 
duces a  random  variable  into  what  might  otherwise  have  been  a 
deterministic  response.  The  on-going  background  activity  masks  and 
distorts  the  ER  such  that  it  usually  is  undetectable  at  the  scalp. 
Thus,  these  two  sources  of  variability  deny  the  investigator  anything 
but  tenuous  conclusions  about  the  relationship  of  a  single  bioelectric 
response  to  specific  stimulus  parameters.  In  order  to  determine  the 
general  characteristics  of  an  ER,  several  data  processing  techniques 
have  been  utilized.  The  most  coinmon  teclinique  is  the  sequential 


average  (see  Qiapter  2)  o£  a  series  of  responses  recorded  from  an 
individual  subject.  Tliis  average  is  fairly  consistent  for  a  parti- 
cular subject  under  identical  stimulus  conditions.  However,  the 
averages  obtained  from  several  subjects  often  appear  to  be  quite 
dissimilar.  Thus,  it  is  difficult  to  establish  a  certain  average 
response  waveform,  elicited  by  specific  stimulus  parameters  across  a 
subject  population. 

In  the  present  study  the  average  of  a  single  response  across  a 
subject  population  is  investigated.  This  average  provides  an  estimate 
of  the  characteristics  of  the  response  to  a  particular  stimulus  which 
are  common  to  an  ensemble  of  subjects.  From  these  common  character- 
istics one  can  develop  a  model  or  template  of  the  average  response 
activity  across  the  ensemble.  The  changes  in  this  response  average 
with  repetitions  of  a  specific  stimulus  can  also  be  investigated.  In 
addition,  the  collection  of  subject  EEG  waveforms  can  be  viewed  as  an 
ensemble  of  independent  sample  functions  from  a  stochastic  process. 
Thus,  assumptions  underlying  the  averaging  process  can  be  examined 
by  classical  tediniques .  These  assumptions  generally  involve  the 
statistical  nature  of  the  evoked  signal  and  background  noise. 

Tlierefore,  the  t^vo  major  objectives  of  this  investigation  are: 

(1)  to  develop  a  better  understanding  of  the  common  cJiaracteristics 
or  norm  of  the  evoked  responses  within  a  subject  population;  and 

(2)  to  provide  an  empirical  technique  to  evaluate  the  prevalent 
statistical  assumptions  about  the  character  of  the  evoked  response 
ajid  background  potentials.  Since  at  this  time  the  response  signal 
cannot  he  separated  from  tlie  background  noise,  the  results  are 


limited  to  some  combination  of  tlie  signal  and  noise  rather  tlian  either 
by  itself.  The  characteristics  of  the  noise  under  non-stimulus  condi- 
tions are  observed,  but  these  may  provide  little  insight  into  the 
background  noise  activity  during  stimulation.  Therefore,  the  second 
objective  is  theoretically  unattainable  at  the  outset.  However,  the 
results  obtained  in  this  study  should  contribute  to  the  knowledge  of 
the  statistical  character  of  evoked  on-going  activity. 

The  results  of  this  investigation  are  divided  into  two  chapters. 
Chapter  3  explores  the  inter- subject  norm  or  ensemble  average  and  the 
transient  clianges  in  this  average.  Also  examined  are  the  non-stimulus 
control  estimates  of  the  average  background  activity.  Chapter  4  con- 
tains an  examination  of  three  assumptions  about  the  signal  and  noise 
which  are  necessary  to  theoretically  justify  results  obtained  by  these 
averaging  techniques.  In  this  examination  the  statistical  character 
of  the  ensemble  of  ER's  with  background  noise  and  the  ensemble  of 
resting  bioelectrical  potentials  is  investigated.  To  facilitate 
understanding  of  the  results  presented  in  these  chapters ,  Chapter  2 
provides  an  overview  of  the  background  material  and  methodology  for 
this  investigation. 


CHAPTER  2 

SEQUEOTIAL  AND  ENSBBLE  AVERAGING 

The  peak-to-peak  amplitude  of  the  potentials  recorded  from  the 
scalp  usually  have  a  range  of  9  to  50  \iv   (Cobb,  1963).  In  most  cases 
the  amplitude  of  the  evoked  activity  in  these  potentials  is  below  the 
threshold  of  observability.  It  therefore  becomes  necessary  for  the 
researcher  to  preprocess  tlie  scalp  potentials  to  extract  some  charac- 
teristics of  the  ER.  Generally,  these  processing  techniques  require 
reiterative  measures  and  estimate  the  response  characteristics  in 
terms  of  probability  statements  or  statistical  averages. 

2.1  Sequential  Averaging 

To  date  the  most  powerful  statistical  descriptor  of  the  ER  is  the 
sequential  average  which  was  implemented  photographically  by  Galambos 
and  Davis  (1943)  and  later  electromechanically  by  Dav;son  (1954)  and 
electronically  by  Barlow  (1957) .  Today  this  average  is  usually  cal- 
culated by  a  general  or  special  purpose  digital  computer.  The  sequen- 
tial average  is  essentially  the  cumulativs.  sum  or  average  of  successive 
trains  of  the  scalp  potentials  which  attend  a  series  of  repeated 
identical  stimuli.  The  mathematical  expression  for  this  average  is 

N-1 
m.(t)=^  E  f.(t+jT)     for  0<t<T  (2.1) 

^     ^^  j=0  ^ 

where  f . (t)  is  the  time  varying  EEG  (electroencephalogram)  recorded 

from  subject  i  and  the  stimulus  occurs  periodically  N  times  at  t=0,  T, 

2T  .  .  . .  The  ~  indicates  the  above  is  an  estimate  of  the  true  mean 


m(t)  for  the  range  of  t  between  0  and  T.  If  the  stimuli  are  visual, 
the  result  of  this  average  is  called  a  visual  evoked  response  (VER) . 
Since  the  present  study  enploys  visual  stimuli,  most  of  the  discussion 
will  focus  on  responses  evoked  in  this  modality. 

The  basic  advantage  of  sequential  averaging  over  other  data  pro- 
cessing m.ethods  is  that  it  extracts  an  estimate  [e.g.,   the  VEK)   of  the 
response  waveshape  from  the  random  potential  fluctuations  monitored 
at  the  scalp.  Other  data  processing  techniques  such  as  autocorrelation 
(Barlow,  1959),  spectral  analysis  (Childers,  et  al. ,  1968),  matched 
filtering  (Negin,  1968) ,  and  synchronous  correlation  (Regan,  1966a) 
may  detect  the  periodicity,  frequency  composition,  presence,  or  phase 
locked  components  of  the  ER,  respectively,  but  they  do  not  produce  the 
waveshape  infoimation  provided  by  averaging.  In  addition,  the  apparent 
underlying  integration  mechanisms  noted  in  most  complex  neurological 
systems  tend  to  justify  or  lend  a  degree  of  rationale  to  the  averaging 
technique,  as  Rosenblith  (1959,  p.  540)  reasons: 

In  the  handling  of  sensory  information,  organisms 
behave  as  if  they  were  acting  on  tlie  basis  of 
activity  averaged  over  substantial  regions  of  the 
nervous  system,  such  a  weighted  averaging  process 
could  be  carried  out  by  comparing  the  outputs  of 
a  large  number  of  neural  elements .  Since  there  are 
at  present  no  devices  of  performing  this  task 
(even  if  one  knew  exactly  where  and  how  to  perform 
it) ,  one  chooses  to  average  over  an  ensemble  of 
responses  to  repeated  identical  stimuli  in  order 
to  bring  out  certain  typical  aspects  of  behavior 
of  the  nei^vous  system. 

The  major  disadvantages  of  sequential  averaging  are:  (1)  only 

those  components  of  the  ER  which  are  time  locked  to  the  stimulus  will 

be  extracted;  (2)  tlie  responses  to  succeeding  stimuli  may  differ  and 

distort  the  average  response  estimate  because  of  the  effects  of  those 


preceding  them;  and  (3)  the  resultant  average  of  maiiy  responses  is 
not  identical  under  invariant  stimulus  conditions.  By  extracting  only 
the  time  locked  components,  the  averaging  process  obscures  any  phase 
and  frequency  deviations  of  the  individual  responses.   In  addition, 
the  degree  of  amplitude  variation  of  these  time  locked  components  is 
hidden  within  the  resultant  waveform.  However,  it  is  impossible  to 
ijTTprove  the  averaging  process  by  weighting  and  phase  locking  each  res- 
ponse without  the  a  priori  knowledge  of  the  specific  angle  or  amplitude 
modulation  which  may  affect  a  series  of  responses.  Recent  studies 
by  Harris  and  Woody  (1969)  and  Woody  (1967)  have  attempted  to  estimate 
the  degree  of  angle  modulation  of  the  ER  by  a  correlation  technique. 
From  this  estimate,  they  have  suggested  an  adaptive  filter  to  improve 
signal  extraction  by  averaging.  However,  to  most  investigators  such 
modulation  is  a  perturbing  factor  distorting  the  estimate  obtained  by 
averaging. 

The  second  disadvantage  is  imposed  by  the  practical  limitation 
that  the  large  number  of  stimuli  necessary  to  extract  an  average  res- 
ponse is  generally  presented  over  a  relatively  short  period  of  time. 
This  limitation  introduces  two  disturbing  factors.  The  first  factor 
can  be  termed  overlap  interference,  i.e.,  the  transient  tail  of  one 
response  overlaps  and  distorts  the  initial  components  of  the  subse- 
quent response.  This  usually  occurs  only  with  relatively  high  stimu- 
lus rates.  The  critical  stimulus  frequency  at  which  this  interaction 
occurs  is  not  clearly  established,  but  the  concensus  of  evidence 
indicates  that  it  lies  in  the  region  of  0.7  Hz  (Barlow,  1960)  to 
2.7  Hz  (Kitzsato,  1966).  Tnere  is  some  indication  that  higher  rates 


are  permissible  if  fjie  stimulus  is  presented  at  random  intervals 
(Ruclikin,  1965). 

The  second  disturbing  factor  is  habituation,  which  can  generally 
be  described  as  the  decrease  in  the  subject's  attention  to  monotonous 
repetitive  stimulation,  which  changes  the  successive  responses.  Since 
the  extraction  by  averaging  of  a  reliable  ER  requires  many  repetitions 
of  tlie  stimulus  (see  Periy  and  Qiilders,  1969)  most  investigators 
discuss  habituation  in  terms  of  the  changes  in  consecutive  averages  of 
large  groups  of  successive  responses.  Therefore  sliort  term  habituation, 
particularly  transient  changes  of  the  ER's  to  the  initial  stimuli,  is 
difficult  to  record.  However,  Zerlin  and  Davis  (1967)  utilizing  a 
highly  responsive  subject,  have  noted  a  systematic  downward  drift  in 
the  amplitude  of  individual  responses  to  successive  auditory  stimuli. 
This  result  agrees  with  the  less  specific  gradual  waxing  and  waaiing 
decrease  in  ajTiplitude  of  consecutive  averages  (\^R's)  found  by  many 
authors  (e.g..  Perry  and  Copenhaver,  1965;  Haider,  et  al. ,  1964;  and 
Garcia-Austt,  et  al. ,  1963).  Bogacz,  et  al.  (1960)  have  also  noted 
that  habituation  to  visual  stim.uli  is  also  reflected  in  increased 
phase  lag  (latency)  of  certain  averaged  response  components. 

2.2  Variability  of  the  VER 

The  third  shortcoming  of  sequential  averaging,  the  variability  of 
tlie  resultant,  poses  the  greatest  dilemma  to  the  researcher  using  this 
tool.  Tlie  XTBR's  taJcen  on  different  occasions  from  a  single  subject 
are  relatively  stable  over  a  period  of  time  (Dustman  and  Beck,  1965); 
however,  when  one  compares  tlie  \'ER's  obtained  under  identical  stim.ulus 
conditions  from  a  population  of  subjects,  variability  of  the  waveform 


makes  it  difficult  to  interpret  the  effect  of  certain  stimulus  para- 
meters on  the  hum^an  visual  system.  Some  investigators  have  attempted 
to  reduce  this  inter-subject  variability  by  expanding  the  resultant 
wavefom  into  orthogonal  sets  for  comparison,  with  little  success 
(see  Donchin,  1966  and  John,  et  al. ,  1964) .  Others  have  organized  a 
system  of  classification  of  the  certain  maxima  and  minima  of  the  VER 
waveform  in  tems  of  amplitudes  and  time  delays  (latencies)  from  stim- 
ulus onset  (see  Bergamini  and  Bergamasco,  1967  and  Kooi  and  Bagchi, 
1964).  With  this  latter  teclinique,  investigators  have  been  able  to 
occasionally  detect  trends  of  diange  common  to  a  subject  population 
with  gross  variations  in  the  stimulus  parameters.  But  such  findings 
are  again  hampered  severely  by  the  variability  of  the  VER.  This 
variability  appears  to  be  the  product  of  changes  in  the  psychological 
state  of  the  subject  (e.g.,  attention)  and  the  inherent  statistical 
character  of  the  BR  and  background  activity  O'-erre  and  Smith,  1964). 

In  an  attempt  to  understand  the  statistical  properties  of  the  ER 
and  background  noise,  a  vast  amount  of  effort  has  been  concentrated 
on  the  classification  of  the  physiological  origin  of  these  scalp 
potentials.  However,  to  date  there  is  only  a  rudimentary  knowledge 
of  the  mechanisms  which  underlie  the  generation  of  these  EEG  wave- 
forms (Elul,  1967).  Tlie  present  evidence  indicates  that  EEG  poten- 
tials arise  from  some  complex  integration  of  excitatoiy  and  inhibitory 
post  synaptic  potentials,  psp  (Humphrey,  1968;  Creutzfeldt  and  Kuhnt, 
1967;  Eccles,  1966  and  Amass i an,  et  al.,  1964).  Tlie  relationship 
between  these  subcranial  potentials  and  the  bioelectric  waveform 
monitored  at  the  scalp  has  been  investigated  .in  animals  by  De  J-ucchi, 


et  al.  (1962) .  They  have  hypothesized  that  the  scalp  acts  as  a 
further  averager  of  psp's.  In  siTidlar  investigations  involving  human 
subjects.  Cooper,  et  al,  (1965)  demonstrate  that  only  an  average  of 
those  potentials  which  are  widely  s>Tichronized  can  be  observed  in  the 
scalp.   It  is  therefore  interesting,  as  Rosenblith  (1959)  notes,  that 
even  with  this  gross  averaging  of  psp's  by  the  scalp  electrodes,  we 
still  find  variability  of  the  ER  to  identical  stimuli. 

Tne  amount  of  tliis  variability  produced  by  the  interaction  of  the 
response  signal  with  the  background  activity  (noise)  is  generally 
unknown.  To  date,  most  investigations  of  this  interaction  have  con- 
centrated on  only  one  conponent  of  the  noise  (resting  EEG) ,  the 
alpha  rhythjTi.  This  component  is  defined  by  Brazier,  et  al.  (1961, 
p.  647)  as  the  "rhythm,  usually  with  frequenc)'  of  8-13  c/sec,  of 
almost  sinusoidal  form,  in  the  posterior  areas,  present  during  relax- 
ation when  the  eyes  are  closed,  attenuated  during  attention,  particu- 
larly visual."  The  physiological  origin  of  this  regular  rhythmic 
activity  is  yet  undetermined.  However,  Andersen  and  Andersson  (1968) 
point  to  the  thalamus  as  the  source  of  the  rhythmic  impulses  which 
trigger  the  psp's  monitored  at  the  scalp.  This  agrees  with  the 
observations  of  Creutzfeldt,  et  al.  (1966).  The  alpha  rhythm  is  the 
m.ost  striking  component  of  the  resting  EEG  of  the  majority  of  normal 
adults.  Tliis  component  is  so  prevalent  that  a  resonant  peak  aroLnid 
10  Hz  can  readily  be  seen  in  the  EEG  frequenc>'  specti^am  O^'^-lter, 
et  al. ,  1967) .  In  addition,  the  late  rhythmic  component  (after 
discharge)  of  the  ER  and  \tR  has  the  general  form  and  dimensions  of 
the  spontaneous  alpha  wave  (Bishop  and  Clare,  1953  and  Barlow,  1960). 
It  is  tiierefcre  not  unexpected  that  studies  involving  tlie  relationship 
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o£  tliis  predominant  EEG  component  to  evoked  responses,  particularly 
tliose  elicited  by  visual  stimuli,  received  precedence. 

A  review  of  the  studies  focusing  on  tJiis  interaction  witli  the  VER 
is  provided  in  the  recent  work  by  Qiilders  and  Periy  (1969) .  They 
note  that  many  investigators  suspect  two  forms  of  8  to  13  Hz. 
activity:  the  first,  the  predominant  non-provoked  EEG  activity,  and 
the  second,  the  components  of  the  evoked  response  Avhich  are  alpha- 
like in  character.  The  evidence  presented  indicates  that  if  two  dis- 
tinct foi-ms  of  activity  exist,  they  probably  share  common  neural 
ensembles.  However,  this  relationship  is  unclear  and  Childers  and 
Perry  stress  the  need  for  further  experimentation  to  secure  the 
additional  information  necessary  to  determine  the  interaction  of 
these  components  and  other  background  activity  with  the  ER. 

2.3  Tlieoretical  Models 

Several  investigators  (e.g.,  Ruchkin,  1969;  Kitasato  and  Hatsuda, 
1965;  Bendat,  1964;  and  Goldstein  and  Weiss,  1962)  attempt  to  under- 
stand the  mechanisms  of  the  VER  and  hence  its  variability  by  developing 
theoretical  models  of  the  statistical  character  of  scalp  potentials. 
These  studies  generally  rely  on  the  h>'pothesis  that  the  signal  (ER) 
and  noise  (background  activity)  are  additive.  This  assumption  pemits 
a  certain  amount  of  mathematical  tractability  to  these  models  but 
contains  little  physiological  basis.  However,  it  is  justified  in 
part  because  of  the  ability  to  extract  a  somewhat  reproducible  esti- 
mate of  the  signal  from  the  scalp  potentials  by  sequential  averaging. 
For  example,  if  the  signal  and  noise  were  multiplicative,  this 


11 


extraction  tedinique  would  provide  conpletely  inconsistent  results. 

In  order  to  facilitate  statistical  analysis  of  the  scalp  poten- 
tials, these  model  developers  generate  an  artificial  stochastic 
process.  This  process  is  constructed  by  segmenting  (usually  at  the 
points  of  stimulation)  the  provoked  EEG  record  into  an  ensemble  of 
shorter  records.  The  statistical  moments  of  this  ensemble  can  then 
be  calculated  from  the  empirical  data  and  compared  to  tiie  theoretical 
model  (e.g.,  the  first  moment  is  the  sequential  average).  Unfortu- 
nately, few  authors  have  progressed  past  the  modeling  stage  to 
compare  the  tlieor>^  with  practice.  This  gap  exists,  because,  as 
Barlow  (1967)  notes,  the  theoretical  ability  to  separate  the  statis- 
tical conponents  (i.e.,  signal  and  noise)  does  not  have  an  experimen- 
tal counterpart.  Such  a  separation  would  require  a  much  greater 
knowledge  of  the  character  of  either  the  signal  and/or  noise  and 
their  interdependence  tlian  is  available  at  the  present. 

In  addition,  the  assunptions  made  about  the  stochastic  model  are 
generally  violated  in  the  empirical  ensemble.  For  example,  the 
empirical  statistics  are  derived  from  temporally  related  members, 
which  requires  a  stationary  ergodic  random  process.  This  is  highly 
improbable  if  one  notes  the  existence  of  a  time  varying  mean  (the 
VER  signal) .  Also  the  members  would  probably  have  some  statistical 
dependence,  reducing  the  quality  of  the  empirical  statistics. 
Finally,  each  subject  would  have  to  be  examined  separately,  which 
would  eliminate  the  estimation  of  the  inter-subject  statistics. 
However,  is  there  a  better  stodiastic  model  and  is  such   a  model 
feasible?  Tliis  question  initiated  tiie  present  study  and  led  to  the 
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development  of  the  following  model  of  a  multi-subject  random  process. 

2,4  N^ilti -Subject  Ensemble  Model 

It  was  decided  that  a  more  appropriate  empirical  stochastic 
process  might  be  generated  from  an  ensemble  of  individual  EEC's 
(i.e.,  each  EEG  which  contains  evoked  activity  is  recorded  from  a 
different  subject).  Such  an  ensemble  would  eliminate  the  temporal 
and  probably  the  statistical  dependence  of  its  members.  Thus,  this 
ensemble  Avould  more  accurately  conply  with  the  classical  definition 
of  a  stochastic  process  as  a  process  composed  of  an  ensemble  of  inde- 
pendent simultaneously  monitored  sample  functions.  The  empirical 
ensemble  statistics  would  then  provide  a  more  precise  evaluation  of 
the  statistical  assumptions  necessar>^  to  theoretically  justify 
averaging  to  extract  a  consistent  response  estimate.  In  addition, 
these  ensemble  statistics,  particularly  the  first  moment  (average) 
of  the  ensemble,  would  pro\ride  information  about  common  inter-subject 
characteristics  of  the  ER. 

The  inlierent  limitations  of  this  artificial  process  are  two-fold. 
First,  the  individual  records  are  not  simultaneously  monitored,  but 
serially  recorded  over  a  several  month  period.  Second,  the  evoked 
acti\dty  in  the  member  functions  may  result  from  several  stochastic 
processes  instead  of  simply  one.  Tlie  first  restriction  becomes 
significant  if  the  nature  of  a  normal  adult  EEG  is  dependent  upon 
tiie  day  of  the  year  or  the  hour  of  the  day.  Fortunately,  however, 
the  present  evidence  points  to  a  relatively  stable  individual  EEG, 
particularly  in  frequency  composition  (Berldiout  and  Walter,  1968  and 
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Jolmson  and  Ulett,  1959).  The  second  limitation  becomes  important  if 
the  population  of  ER's  can  be  divided  into  discrete  classes.  To 
date,  there  has  been  some  indication  (Perry,  et  al. ,  1968)  that  such 
sub-divisions  might  exist;  however,  the  method  of  classification  and 
the  number  of  classes  is  yet  undetermined.  Therefore,  the  present 
research  effort  will  assume  a  continuum  of  EEG  and  response  types 
which  can  be  modeled  under  the  same  stochastic  process. 

2.5  Ensemble  Averaging 

The  ensemble  mean  or  average  provides  an  estimate  of  the  common 
inter-subject  response  to  a  single  stimulus.  The  mathematical 
expression  for  this  average  is 

1  N 
m.Ct)  =  i  Z  f  (t)  (2.2) 

^     ^^  i=l  ^ 
where  f - (t)  is  the  EEG  containing  evoked  activity  recorded  from 
subject  i  out  of  a  population  of  N  subjects.  This  ensemble  average, 
VER,  and  unprocessed  scalp  potentials  are  schematically  represented 
in  Figure  2.1.  This  figure  demonstrates  the  basic  difference  in  the 
ensemble  and  sequential  average  (i.e.,  \T:R) .  The  ensemble  average 
provides  an  estimate  of  the  inter-subject  average  response  to  a  single 
stimulus  and  the  \T:R  provides  the  inter- stimulus  average  response  for 
a  single  subject.  The  ensemble  average  eliminates  the  inter-subject 
response  variability  while  it  reflects  the  changes  in  the  average 
response  to  a  series  of  individual  stimuli.  Ihe  WR   average  elimi- 
nates the  inter- stimulus  response  variability  while  it  reflects  the 
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individuality  of  a  subject's  average  response. 

Hie  average  of  a  set  of  \'ER's  from  a  subject  ensemble  equals  the 
average  of  the  set  of  ensemble  averages  for  a  series  of  stimuli,  pro- 
vided both  averages  are  obtained  from  the  same  inclusive  collection 
of  data.  This  double  average  provides  an  overall  estimate  of  the 
response  activity  by  sacrificing  all  variability  infoimation.  This 
overall  average  has  been  suggested  by  Da\\'Son,  et  al.  (1968)  to  reduce 
all  excess  infoimation  about  the  response.  Therefore,  this  single 
estimate  might  provide  a  general  inter-subject,  inter- stimulus 
response  norm  to  specific  stimulus  parameters. 

2.6  Methods 

In  the  initial  stages  of  this  investigation  it  was  decided  that 
an  ensemble  size  (subject  population)  of  fifty  would  provide 
reasonable  statistical  accuracy  and  remain  within  the  realm  of  experi- 
mental practicability.  Both  stimulus  (signal)  and  non-stimulus 
(noise)  trials  were  recorded  from  each  of  the  subjects.  The  stimulus 
was  a  monoailar  periodic  flash  of  light  which  had  a  repetition  rate 
of  about  one  flash  per  second.  The  stimulus  was  visually  occluded 
during  noise  (control)  trials.  The  potentials  recorded  in  these 
trials  were  monitored  from  a  pair  of  scalp  electrodes  vertically 
oriented  on  the  occiput  midline.  These  potentials  were  amplified 
by  a  factor  of  100  K;  separated  into  t^vo  data  diannels  (i.e.,  a 
wide  band,  0.3  to  50  Hz,  channel  and  a  narrow  band,  0.3  to  15  Hz, 
channel);  and  recorded  on  magnetic  tape. 

These  recorded  potentials  were  then  digitized  and  stored  in  an 
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IBM  1800  computer  in  order  to  facilitate  necessary  mathematical 
computation.  Computer  memory  limitations  restricted  the  amount  of 
data  analyzed.  Therefore,  most  of  the  results  presented  in  this 
study  are  derived  from  activity  following  ten  stimulus  triggers  (both 
signal  and  noise  trials)  selected  from  the  periodic  (approximately 
1  Hz)  stimulus  train.  The  triggers  chosen  are  the  first  through 
fifth,  tenth,  twentietli,  thirtieth,  fortieth  and  fiftieth  in  the 
repetitive  stimulus  train.  All  of  the  computation  was  performed  on 
the   first  950  milliseconds  of  data  following  each  trigger.  A  more 
detailed  discussion  of  the  data  collection  and  analysis  procedure 
is  provided  in  i^pendix  A. 

2.7  Objectives 

Tlie  two  objectives  of  this  study,  as  stated  in  the  introductory 
chapter,  are:  (1)  to  develop  a  better  understanding  of  the  common 
characteristics  or  norm  of  the  evoked  responses  within  a  subject 
population  and  (2)  to  provide  an  empirical  technique  to  evaluate 
the  prevalent  statistical  assumptions  about  the  diaracter  of  the  ER 
and  background  potentials.  Tlie  inter-subject  ensembles  recorded  under 
stimulus  and  non-stimulus  conditions  provide  data  from  which  both 
these  objectives  can  be  pursued.  The  average  across  the  ensemble 
under  stimulus  conditions  furnishes  an  estimate  of  the  common  inter- 
subject  evoked  activity.  The  data  provided  by  these  ensembles  is 
appropriate  for  the  investigation  of  theoretical  assumptions  about 
the  EEC  stochastic  process  since  these  potentials  are  a  reasonable 
facsimile  of  such  a  process. 
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In  pursuit  of  the  first  objective  the  ensemble  average  across 
fifty  subjects  under  stimulus  conditions  is  examined.  It  is  antici- 
pated that  the  results  from  this  average  will  reveal  the  components  of 
the  evoked  activity  which  are  the  most  consistent  across  the  subject 
population.  These  components  should  reappear  after  each  stimulus  with 
some  degree  of  variability.  The  variability  would  probably  be  the 
result  of  transient  changes  (e.g.,  habituation)  in  tlie  inter-subject 
response  with  the  nimiber  of  stimuli,  variations  introduced  by  the 
inherent  error  of  the  estimation  process,  or  background  noise  poten- 
tials which  may  remain  in  the  average.  A  transient  variation  should 
result  in  some  overall  trend  of  change  in  the  response  average  while 
noise  interference  and  estimation  error  should  affect  the  average  in 
a  more  random  manner.  The  transient  trends  of  response  estimate  and 
confidence  in  this  estimate  across  a  series  of  stimuli  will  be  inves- 
tigated. Tlie  ensemble  average  under  non-stimulus  conditions  is 
expected  to  provide  an  approximation  of  the  amount  of  noise  contained 
in  the  response  estimate.  The  average  of  a  collection  of  ensemble 
averages  should  eliminate  transient  variability,  reduce  the  components 
of  the  backgroand  noise  and  provide  a  fairly  stable  inter-subject 
response  template. 

In  order  to  theoretically  justify  the  above  response  estimation 
by  averaging,  it  is  necessary  that  (1)  the  background  noise  come  from 
a  relatively  wide  sense  stationary  random  process ;  (2)  the  signal  and 
noise  processes  be  statistically  independent;  and  (3)  the  signal  and 
noise  be  additive.  Since  the  signal  and  noise  potentials  cannot  be 
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separated  at  the  present  time,  the  above  properties  are  generally 
assumed  to  be  true.  The  second  objective  of  this  study,  evaluation 
of  these  assumptions,  is  hindered  by  the  inability  to  segregate  the 
signal  and  noise.  However,  by  hypothesizing  one  or  two  of  these 
assumptions  an  attempt  can  be  made  to  empirically  assess  the  third. 
For  example,  by  hypothesizing  the  signal  and  noise  to  be  statistically 
independent,  the  inter-subject  ensemble  under  non-stimulus  conditions 
should  approximate  the  noise  activity  under  stimulus  conditions.  The 
stationarity  (i.e.,  time  dependence  of  the  ensemble  statistics)  of 
the  background  noise  can  then  be  examined  using  the  resting  EEG  as 
noise  approximation.  Evaluation  of  the  second  assumption  can  be 
obtained  in  a  similar  manner  by  hypothesizing  the  third  assumption 
and  examining  the  ^^  ensemble  average  and  ensemble  variance  under 
stimulus  and  non-stimulus  conditions.  The  assessment  of  the  third 
assumption  is  more  difficult.  However,  by  hypothesizing  the  first 
two  assumptions  we  can  examine  tiie  growth  of  the  signal  to  noise 
CS/N)  ratio,  speculating  the  signal  to  be  the  ensemble  average  and 
tlie  noise  the  +_  ensemble  average  under  stimulus  conditions.  If  the 
signal  and  noise  are  additive,  the   S/N  ratio  should  grow  linearly 
with  the  nmnber  of  averages. 

The  results  presented  in  the  two  succeeding  chapters  attempt  to 
realize  the  above  objectives.  Tlie  response  activity  common  to  tlie 
subject  ensemble  is  first  examined,  followed  by  the  evaluation  of  the 
statistical  assuinptions  about  the  diaracter  of  the  ER  and  background 
noise. 


QiAPTER  3 
RESULTS  OF  ENSBBLE  AVERAGING 

3.1  Introduction 

In  this  chapter  the  inter-subject  average  response  to  a  monocular 
visual  stimulus  is  examined.  The  average  across  the  fifty  subject 
ensemble  record  (see  Figure  2.1)  is  utilized  to  extract  an  estimate 
of  the  response  activity  from  tlie  random  scalp  potentials .  This 
ensemble  average  should  provide  an  inter-subject  response  norm  if: 
(1)  the  assumptions  mentioned  in  the  previous  chapter  hold;  (2)  there 
exist  stimulus  locked  response  components  which  are  common  to  the 
subject  population;  and  (3)  tlie  amplitude  of  the  average  background 
activity  in  this  estimate  is  much  less  tlian  the  amplitude  of  these 
response  components.  The  first  condition,  tlie  validity  of  assumptions 
about  the  response  signal  and  background  noise  necessary  to  justify 
averaging,  is  the  subject  of  Chapter  4.  The  present  chapter  will 
hypothesize  these  assumptions  to  be  true  and  examine  the  ensemble 
average  for  any  inter-subject  stimulus  locked  response  components.  In 
addition,  an  estimate  of  the  amount  of  background  noise  contained  in 
this  average  is  investigated. 

In  order  to  simplify  discussion  of  the  ensemble  average  evoked 
response  estimate  the  term  ensemble  evoked  response  (EER)  will  be 
utilized.  Ihis  term  is  analogous  to  tlie  term  visual  evoked  response 
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(\'ER)  in  that  it  is  not  intended  to  iiiply  that  the  average  response 
estimate  contains  only  evoked  activity.  Theoretically  the  \T:R  and  EER 
contain  noise  coinponents  whose  variance  decreases  directly  with  the 
number  o£  sums  in  the  average  (see  Perry  and  Qiilders ,  1969) .  An 
estimate  of  the  amount  of  noise  in  the  EER  is  obtained  from  the  average 
across  the  ensemble  under  non-stimulus  conditions.  Tliis  average  pro- 
vides a  series  of  control  waveforms,  each  attending  a  particular 
occluded  stimuli,  whose  amplitudes  should  provide  an  estimate  of  the 
background  noise  in  the  EER,  This  background  noise  activity  which  is 
not  related  to  the  stimulus  will  introduce  variations  in  the  EER  from 
one  stimulus  to  the  next. 

As  noted  in  the  previous  chapter,  other  parameters  which  could 
affect  the  EER  across  stimuli  are  doAaations  of  the  estimate  about 
tlie  true  mean  and  variations  in  this  response  mean  with  stimulus 
number.  Tlie  deviation  of  the  EER  from  the  true  mean  is  tlie  result  of 
the  variance  of  the  potentials  across  the  ensemble  (see  Bendat  and 
Piersol,  1966).  Both  the  response  signal  and  background  noise  contri- 
bute to  this  variance  since  it  is  unlikely  that  the  response  activity 
from  every  subject  is  identical  and  the  noise  is  random.  However, 
since  the  signal  activity  generally  has  less  amplitude  (Perry  and 
Childers,  1969)  and  greater  similarity  across  the  ensemble  than  the 
noise  potential,  the  latter  should  account  for  most  of  the  ensemble 
variance.  Therefore  the  estimation  deviations  of  the  EER  should  be 
primarily  the  result  of  tlie  background  noise.  The  range  of  these 
deviations  can  be  predicted  utilizing  the  Tchebycheff  inequality  and 
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ensemble  variance.  Tliis  prediction  takes  tlie  form  of  confidence 
intervals  about  the  EER. 

In  order  to  examine  the  second  parameter,  the  variations  in  the 
EER  caused  by  changes  in  the  actual  inter- subject  response,  it  is 
necessary  to  separate  response  signal  from  the  noise  activity  in  the 
ensemble  average.  Since  this  is  not  possible  at  this  time,  the 
present  investigation  will  generally  be  restricted  to  variations  in 
the  more  prominent  components  of  the  EER's  which  are  fairly  consistent 
across  stimuli.  Since  the  stimuli  are  identical,  these  variations 
should  result  in  some  trend  of  change  rather  than  random  fluctuations. 


3.2  Selection  of  Bandt'/idth 

The  EER  and  other  results  presented  in  this  chapter  and  the 
succeeding  chapter  are  obtained  from  the  stimulus  provoked  and  resting 
EEG  ensembles  mentioned  in  Chapter  2.  In  addition,  most  of  the  dis- 
cussion of  these  results  is  limited  to  narrow  band  filtered  (0.3  to 

Although  the  fidelity  of  this  more  restrictive 
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Bandwidtli  Comparison 

Tiie   ensemble  mean  ajid  standard  deviation 
for  tis'O  bandwidths  are  superimposed. 
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bandwidth  is  not  as  good  as  tiie  wider  band  potentials  (0.3  to  50  Hz), 
it  still  retains  the  significant  EHG  potential  fluctuations,  as 
Figure  3.1  demonstrates.   In  this  figure  the  first  and  second  ensemble 
moments  for  both  bandwidths  are  superimposed  and  reflect  only  minor 
variations,  including  a  slight  phase  shift  introduced  by  the  filter. 
In  addition,  the  upper  limit  of  15  Hz  exceeds  the  minimum  range  of 
ten  times  the  frequency  of  stimulation  suggested  by  Qiilders,  et  al. 
(1968)  to  recover  a  \'TSR  average  witli  good  fidelity.  Therefore,  only 
the  narrow  band  averages  which  provide  an  accurate  description  of  the 
EEG  activity  will  usually  be  discussed.  However,  many  of  the  results 
presented  here  are  duplicated  for  the  wider  band  potentials  in 
Appendix  B,  This  appendix  also  contains  other  results  which  are  not 
pertinent  to  the  discussion.  Any  reference  to  these  results  will 
have  a  prefix  B. 

3.3  Ensemble  Evoked  Response 

The  ensemble  evoked  response,  EER,  to  the  selected  photic  stimuli 
is  provided  in  Figure  3.2.  Each  average  is  designated  by  the  number 
signifying  which  stimulus  within  the  train  evoked  it.   In  each 
waveform  the  flash  stimulios  occurs  at  the  extreme  left  of  each  950 
millisecond  response  average.  Ihe  first  four  ensemble  response 
averages  can  also  be  found  in  Figure  2.1,  which  depicts  the  method  by 
which  the  averages  were  obtained.  Tlie  wavefonns  are  segmented  into 
seven  135  millisecond  intervals  to  facilitate  comparison  and  dis- 
cussion. Again,  it  should  be  pointed  out  that  eacli  EER  is  only  an 
estimate  of  the  actual  evoked  potential  and  contains  both  evoked  and 
background  activity. 
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Interval 
1  2  3  4  5  6  7 


FIGURE  3.2 

Ensemble  Evoked  Response 

Ensemble  averages  for  ten 
selected  stimuli  are 
divided  into  7  intervals 
or  135  irdlliseconds . 


A  cursory  glance  at  the  EER's  in 
Figure  3.2  provides  the  reader  with  little 
indication  that  the  fluctuating  waveforms 
contain  any  response  activity.  Havever, 
if  one  utilizes  the  inverted  sinusoid  in 
the  second  interval  of  the  fourth  res- 
ponse as  an  ER  template,  the  most  common 
attribute  of  the  wavefoim  set  becomes 
more  apparent.  This  inverted  sinusoid, 
although  sometimes  partially  distorted, 
appears  in  the  second  interval  of  each 
response  in  the  figure.  The  potential 
fluctuations  across  the  ten  EER's  also 
share  a  degree  of  likeness  in  the  latter 
portion  of  the  first  and  initial  part  of 
the  third  intervals .  However,  neither 
approaches  the  conformit)''  found  in  the 
second  interval. 

In  order  to  better  establish  the 
bounds  of  the  response  activity  it  is 
necessar)^  to  find  some  quantitative 
measure  of  this  similarity  across  the 
EER's.  The  degree  of  simlarity  between 
two  responses  for  a  specific  interval 
can  be  quantified  by  a  correlation 
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coefficient  which  would  provide  a  total  of  45  correlations  per 
interval  for  the  ten  EER's  of  Figure  3.2.  This  coefficient,  p,  can 
be  expressed  as: 

D  =  -^^  (3.1) 

1   ^x  y 

where  H'       is  the  time  covariance  between  EER's  x  and  y  in  interval 
xy 

i  and  ^     and  T     are  the  time  standard  deviations  of  these  wavefonns 
X  y 

in  this  interval.     How  to  reduce  these  paired  correlations  to  a  single 
measure  of  the  similarity  across  the  entire  set  of  EER's  becomes 
the  major  concern.     The  method  chosen  in  this  study  was  to  average 
the  correlation  coefficient  for  one  interval.     Tliis  average  provides 
a  conservative  measure  of  the  conformity  across  the  ten  responses 
since  most  of  the  correlation  coefficients  betu^een  waveform  pairs 
must  be  significant  to  obtain  a  significant  average  correlation.     The 
average  correlation  reduces  the  effects  of  sporadic  deviations  in 
these  correlation  coefficients ,  which  randomly  appear  between 
paired  wave foim  segments.     The  significance  of  tiiese  correlations  is 
questionable,  since  they  do  not  have  any  consistent  relation  to  the 

stimulus. 

The  correlation  coefficients  for  each  paired  EER  (Figure  3.2)   at 
a  specific  intewal  are  provided  in  Tables  B.l  through  B.7  in 
Appendix  B.     Tlie  averages  and  maxima  and  minima  of  these  correlations 
are  provided  in  Table  3.1.     Tlie  average  correlations  reinforce  the 
previous  observ^ation  that  the  activity  in  the  second  intei-val  across 
the  ten  EER's  is  the  most  consistent.     In  fact,  even  the  minimum 
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TABLE  3.1 
Average  Correlations  for  Figure  3.2 


1 

2 

3 

4 

5 

6 

7 

30 

.88 

.12 

.05 

.06 

.04 

-.07 

89 

.99 

.88 

.88 

.94 

.87 

.98 

45 

.58 

-.71 

-.92 

-.87 

-.92 

-.87 

INTERVAL 

p 

p 

^min 


correlation  of  .58  is  above  the  1  per  cent  significance  level,  i.e., 
assuming  samples  are  normally  distributed  and  independent  (Underwood, 
et  al. ,  1954) .  The  average  of  the  correlation  coefficients  across 
the  segmented  EER's  obtained  from  the  wider  band  potentials  were 
generally  less  than  those  in  Table  3.1  (see  Tables  B.8  through  B.13). 
Thus  filtering  tlie  data  appears  to  eliminate  some  of  the  variability 
across  the  response  estimates. 

The  first  interval  of  the  EER's  is  the  only  other  135  millisecond 
time  period  where  the  waveforms  demonstrate  any  degree  of  consistent 
correlation.  This  consistency  appears  to  be  the  result  of  the  more 
common  components  in  the  latter  half  of  this  interval  rather  than 
the  initial  deflections  after  stimulus  onset.  In  intervals  3  through 
7  (i.e.,  270  to  950  milliseconds  after  the  stimulus)  the  EER's 
demonstrate  negligible  correlation.  This  variability  of  tliese  late 
conponents  of  the  EER  seems  to  indicate  that:  (1)  the  average  back- 
ground activity  is  much  greater  than  tlie  average  response  activity 
during  this  period;  (2)  the  individuality  of  the  response  for  different 
subjects  is  greater,  whicli  increases  tlie  deviations  of  tlie  estimate 
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in  tliis  region;  or  (3)  the  variations  with  stimulus  number  are  greater 
in  tlie  late  response.  In  similar  correlation  tests,  Perry  and 
Childers  (1969)  have  demonstrated  a  significant  decrease  in  the 
average  correlation  coefficients  of  inter- subject  VER's  following 
the  first  300  milliseconds  of  the  response,  and  a  decrease  in  intra- 
subject  correlations  after  500  lailliseconds .  They  also  note  that 
the  highest  correlations  are  generally  found  in  the  second  and  third 
100  millisecond  intervals,  which  is  in  close  agreement  with  the 
previous  results.  Tlierefore,  the  EER  and  VER  results  seem  to  reject 
the  reasoning  that  it  is  the  individuality  of  the  evoked  response 
which  introduces  the  variability  in  the  late  components  of  the 
ensemble  average. 

It  is  interesting  to  note  the  occurrence  of  alpha-like  activity 
in  the  latter  portion  of  response  40.  This  singular  occurrence  may 
be  the  result  of  a  particular  combination  of  the  random  EEG  potentials. 
However,  these  potential  fluctuations  have  a  strong  resemblance  to 
tlie  after- discharge  activity  noted  in  the  studies  mentioned  in 
Chapter  2. 

In  Figure  3.3  the  ten  EER  waveforms  are  superimposed  with  the 
time  mean  removed  from  each  average.  In  this  figure  the  consistency 
of  the  responses  in  the  region  from  about  60  to  300  milliseconds 
after  the  stimulus  can  easily  be  seen.  Because  of  the  similarity  of 
tlie  EER  in  this  region  (R^)  it  will  be  designated  the  primary 
response  region.  In  this  region  there  are  three  prominent  peaks. 
These  peaks  are  identified  in  Figure  3.3  as  components  one,  two 
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aiid  three,  and  occur  approximately  90  msec,  150  msec,  and  250  msec 
after  stimulus  onset,  respectively. 

The  average  correlation  coefficient  in  the  primary  response 
region  is  0.75,  whidi  is  slightly  less  than  the  average  correlation 
in  the  second  interval  of  Figure  3,2.  This  decrease  is  probably  the 
result  of  the  increase  in  interval  width  from  135  to  240  msec.  Thus, 
the  second  region  of  Figure  3.3  has  more  components  and  greater 
variability  than  the  second  interval  of  Figure  3.2.  The  average 
correlation  of  the  tliird  region  in  Figure  3.3  is  0.23,  which  indi- 
cates a  minor  degree  of  similarity  not  previously  found  in  the  corres- 
ponding correlations  in  Table  3.1.  This  similarity  is  possibly  the 
result  of  some  low  frequency  activity  common  to  the  tails  of  the  EER's 
which  would  not  be  detected  by  the  correlations  of  the  135  millisecond 


FIGURE  3.3 

Superimposed  Ensemble  Evoked  Response 

Overlapping  EER's  are  divided  into  three  regions  for  comparison. 
R^  =  0  to  60  msec;  Y^^   =  60  to  300  m^ec;  R^  =  300  to  950  msec. 
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intei'vals.  However  such  a  low  frequency  component  is  not  visually 
detectable,  except  for  a  very  slight  negative  drift  in  tJie  late 
response  components. 

3.4  Noise  in  the  Enseinble  Evoked  Response 

The  EER's  in  the  previoiis  figures  contain  an  unknown  amount  o£ 
background  noise.  One  estimate  of  the  amplitude  of  the  background 
activity  vsrhich  remains  in  the  EER  is  the  ensemble  average  under 
non-stimulus  (control)  conditions.  The  accuracy  of  this  estimate  is 
dependent  on  the  relationship  between  the  signal  and  the  noise.  For 
example,  if  tlie  signal  and  noise  are  independent  and  the  noise  is 
stationary,  this  average  provides  a  good  measure  of  the  amount  of 
background  activity  in  the  EER.  The  waveshape  of  this  average  is 
important  only  if  a  particular  component  appears  time  locked  to  the 
occluded  stimulus,  indicating  the  existence  of  some  experimental 
artifact . 

In  Figure  3.4,  this  average  is  subdivided  into  ten  control  wave- 
forms attending  the  ten  occluded  stimuli.  The  format  of  this  figure 
is  identical  to  that  of  Figure  3.2.  The  average  correlations  of  the 

TABLE  3.2 
Average  Correlations  for  Figure  3.4 
INTERVAL        12      3      4      5 
p,,,„  .07    -.02    -.03     .01    -.07 

P,„3_^  .81     .92     .95     .99     .87 

Pj^^  -.79    -.99    -.93    -.92    -.98 
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FIGURE  3.4 
Ensemble  Control  Averages 

Ensenil-ile  averages  for  ten 
occluded  stimuli  are 
divided  into  7  ijitervals 
of  135  milliseconds. 


intervals  of  Figure  3.4  given  in  Table  3.2 

demonstrate  the  lack  of  correlation  between 

the  control  waveforms.  These  correlations 

are  similar  to  the  average  EER  correlations 

in  the  last  four  intervals  of  Table  3.1. 

This  supports  the  previous  observation 

that  there  may  be  little  response  activity 

in  this  region. 

In  every  interval  the  root  mean  square 

of  the  control  waveforms ,  C_^^ ,  is  less 

'  rras 

tlian  the  vms   value  of  the  responses,  R   > 

in  the  second  interval  of  Figures  3.2 

(Avhere  these  rms  values  have  the  dc  removed) 

The  minijmum  R    is  only  1.03  times  the 
rms      ^ 

maximum  C    and  the  magnitude  of  the 
rms         ^ 

largest  R    is  9.7  times  the  smallest 
^  rms 

C   .  The  average  of  all  ratios  of  R , 

rms  °  rms 

to  C    is  2.6.  Therefore,  assuming  that 
rms  '       •=" 

the  ensemble  control  average  provides  a 
good  estimate  of  the  amount  of  background 
activity  remaining  in  the  averaged  response 
and  that  the  signal  and  noise  are  additive, 
what  is  the  signal  to  noise  ratio  (S/N)  in 
the  EER?  To  date,  the  answer  to  this 
question  cannot  be  resolved,  since  it  is 
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necessary  to  also  kr.ow  the  phase  o£  the  average  background  activity. 
However,  with  only  a  magnitude  estimate,  the  upper  and  lower  bounds  of 
this  signal  to  noise  ratio  can  be  found  by  adding  and  subtracting  the 
magnitude  of  the  control  from  the  EER  and  dividing  by  the  control 
amplitude,  i.e. , 


^rms  "  ^rms  <  g/N  <  ^™^  ""  ^"^         (3-2) 
r      ~    ~     C 

It  should  be  noted  that  in  the  case  when  the  lower  bound  is  negative 
(which  is  a  theoretically  inpossible  event)  the  lo^s'er  bound  is  set 

to  zero. 

Rather  than  establish  these  bounds,  a  more  common  practice  is 
to  estimate  the  signal  to  noise  ratio  in  tiie  average  response  by 
dividing  tlae  magnitude  of  the  average  response,  R^,  by  the  average 
noise,  C^^  (Perr>',  1966),  i.e., 

S  ^  ^TT1S_  (3.3) 

Using  this  estimate,  the  relative  signal  to  noise  measures  on  the 
previous  page  provide  tiie  maximum,  minimum  and  average  S/N.     Hius,  one 
would  predict  that  the  EER's  in  the  second  interval  of  Figure  3.2 
contain  approximately  10  to  50  per  cent  noise  with  an  expected  value 
of  28  per  cent.      In  Chapter  4  it  is  sho;vn  that  under  certain  assunp- 
tions  these  predicted  S/N  values  are  reinforced  by  the  degree  of 
similarity  (i.e.,   correlation)  previously  fo'ond  bet^veen  tiie  EER's. 
Therefore,  the  average  amplitude  under  control  conditions  probably 
provides  a  reasonable  estimate  of  the  a^Bount  of  noise  distorting  the 
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EER  waveform. 

I£  the  overall  rms  value  of  the  control  waveforms  (approximately 
0.5  yv)  represents  the   expected  noise  distortion  in  the  response 
average,  then  Figure  3.5  provides  an  upper  and  lower  bound  of  the 
noise  distortion  expected  in  the  EER.  The  actual  magnitude  of  the 
background  noise  potentials  at  any  instant  could  be  several  times 
tiieir  rms  value.  However,  the  bounds  in  Figure  3.5  provide  an 
average  approximation  of  the  amplitude  of  the  background  noise  which 
could  be  contained  in  the  EER  estimate.  That  is,  the  average 
response  without  noise  should  be  expected  to  be  within  these  bounds 
if  the  control  waveforms  provide  a  good  estimate  of  background 
noise  and  if  the  EER  is  a  good  estimate  of  the  average  signal  plus 
noise  activit)^.  The  latter  stipulation,  the  quality  of  the  EER 
estimate,  provides  further  confidence  limits  in  the  average  evoked 
activity.  Thus  a  more  accurate  bound  on  the  ei^oked  activity  is  a 
combination  of  the  noise  bounds  and  the  confidence  limits  on  the 
EER. 

^f^^  ^/m^   ^^%  i^H^  /^^ 

12  3  4  5 

Response  Number 
10  20  30  40  50 

^^W^  m^m  u^%4  ^v^  ^1^^ 

FIGURE  3.5 

Noise  Amplitude  in  EER 

Ov^erall  rsTis  control  value  is  added  and  subtracted  from  the  center 
EER  waveform. 
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3.5  Quality  of  the  Ensemble  Evoked  Response  Estimate 

llie  EER  is  a  good  estrmator  o£  the  signal  plus  noise  mean  if  the 
deviations  of  this  estimate  about  the  mean  are  relatively  minor.  A 
measure  of  tJiis  estimation  deviation  is  provided  by  the  standard 
deviation,  a  .(t) ,  of  the  underlying  process.  The  actual  value  of 
OrCt)  is  uhkno^m.  However,  an  unbiased  estimate  of  oAt)   is  the 
ensemble  standard  deviation,  ar(t) ,  which  is  calculated  from  the 
equation 


1   N   .  2 


1/2 


(3.4) 


where  N  is  the  number  of  f. (t)'s,  subject  EEC's,  and  m^Ct)  is  the 
ensemble  average  given  by  equation  2.2. 

A  plot  of  aj-(t)  for  the  potentials  attending  each  of  the  ten 
is  provided  in  Figure  B.5.  These  ten  waveforms  are  superimposed  in 
Figure  3.6  in  order  to  observe  any  stimulus  effects  on  a^Ct) .  As 

5-r 


FIGURE  3.6 
Superimposed  Ensemble  Standard  Deviations 


Each  standard  deviation  attends  one  of  the  ten 
selected  stimuli  which  occur  at  the  extreme 
left  of  the  figure. 
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can  be  seen  in  this  figure,  the  fluctuations  of  the  a^(t)  waveforms 
do  not  appear  to  have  any  consistent  relationship  to  tlie  stimulus, 
other  tlian  a  slight  average  increase  around  the  primary  response 
region. 

If  the  ensemble  standard  deviations  in  Figure  3.6  provide  a  good 
estimate  of  the  standard  deviations  of  the  process,  then  the  standard 
deviation  of  the  mean  estimate,  mr(t) ,  is 

sd[m.(t)]  ^  — (3.5) 

^       v^ 

This  standard  deviation  can  be  incorporated  into  the  Tchebycheff 

inequality  to  provide  reliability  or  confidence  bounds  on  the  EER 

estimate.  This  inequality  provides  a  lo\\'er  bound  on  the  probability 

that  the  true  mean  and  estimated  mean  will  differ  by  less  than  n 

standard  deviations  of  the  estimate.  For  example,  we  can  be  at  least 

89  per  cent  confident  that  the  true  mean  and  estimated  mean  (EER) 

will  differ  by  no  more  than  three  standard  deviations  (calculated 

by  Equation  3.5)  of  the  estimate,  regardless  of  the  underlying 

probability  distribution.  If  the  underlying  probability  distribution 

is  Gaussian,  the  confidence  increases  to  99  per  cent. 

Tlie  89  per  cent  confidence  intervals  for  ten  EER's  are  plotted 

in  Figure  3.7.  The  boimds  in  this  figure  are  constioicted  from  ten 

EER  waveforms  plus  or  minus  three -sevenths  of  the  ensemble  standard 

deviations  superimposed  in  Figure  3.6.  As  can  be  seen  in  Figure 

3.7,  or  calculated  from  Figure  3.6,  the  overall  width  of  the  intervals 

varies  from  about  2  to  4  pv.  Although  this  appears  to  be  a  relatively 
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large  amount  of  variability,  the  shape  o£  tlie  bounding  waveforms 
generally  conforms  to  the  EER  waveshape,  with  a  few  exceptions.     The 
most  notable  exception  occurs  in  the  first  response,  around  300 
milliseconds  after  stimulus  onset.     This  variation  could  have  been 
anticipated  from  Figure  3.6,  since  tlie  standard  deviation  of  the 
potentials  in  this  region  of  the  first  response  partially  stands 
above  the  remaining  response  deviations. 

If  the  confidence  intervals  in  Figure  3.7  are  compared  with  the 
■VER  confidence  intervals  in  Figure  3.8,   the  reliability  of  the  'V'ER 
is  greater  than  the  EER  in  two  out  of  the  three  cases  sho^vn  (J.M. 
and  B.A.).     In  the  third  case   (J.F.)   the  reliabilities  are  approxi- 
mately equal,   if  the  greater  amplitude  of  the  VER  is  taken  into 
consideration.      In  all  three  cases  the  XTER's  were  obtained  by 
averaging  the  subject's  response  to  fifty  stimuli.     Therefore,  for 
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FIGURE  3.7 
EER  Confidence  Intervals 


The  above  intervals  are  plus  and  minus  three  standard  deviatioiis  of 
the  EER  estimate.  This  is  about  3/7  of  the  deviations  in  Figure  3.6. 
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the  same  muriber  of  averages,  the  EER  appears  to  provide  as  good  an 
estimate  of  the  response  activity  as  some  \'ER's. 

The  \^R  amplitude  is  about  one  and  one-half  times  that  of  the 
EER's  in  all  of  the  cases  shavn.  Tliis  larger  amplitude  and  greater 
reliability  of  the  \'ER  estimate  implies  that  the  inter-subject 
response  variability  is  greater  than  the  variability  of  responses 
obtained  from  a  single  subject.  This  is  what  one  would  predict  from 
the  comparison  of  VER's  obtained  from  a  single  subject  and  those 
obtained  from  several  subjects.  In  addition,  the  background  noise 
variability  would  influence  the  reliabilit)^  of  the  estimate.  Thus 
inter-subject  EEG  variability  (e.g.,  amplitude  and  frequency  deviation 
in  the  background  activity)  probably  reduces  the  confidence  in  the 
EER  estimate.  One  variable  parameter  which  might  reduce  confidence 
in  the  \^R  and  not  in  the  EER  would  be  a  general  change  in  the  evoked 
activity  witli  the  number  of  stimuli.  In  order  for  this  variation  not 
to  affect  the  EER  it  is  necessary  for  this  change  to  be  fairly  con- 
sistent across  the  stimuli  for  tlie  subject  population. 
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FIGURE  3.8 

\T:R  Confidence  Inter\''als 

Scale  and  procedui'e  is  the  same  as  in  Figure  3.7.     Each  of  tlie  above 
subjects  are  members  of  the  ensemble. 
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3.6     Vari^ti^n^  in  tiie  Ensen^ 

As  noted  previously,   the  EER  provides  an  estimate  of  the  changes 
in  the  evoked  response   (comon  to  tiie  subject  population)  with  repe- 
titions o£  tiie  stimulus.     Tliat  is,  each  EER  in  Figure  3.2  is  an 
estimate  of  the  activity  which  is  time  locked  to  a  specific  stijnulus 
witliin  the  periodic  train.     In  order  to  obtain  an  accurate  estimate 
of  the  response  changes  related  to  tlie  stimulus  nuiii^er,  it  is  necessary 
to  remove  the  EER  variations  resulting  from  the  background  noise. 
Since  this  cannot  be  accomplished  at  this  time,  the  EER's  across 
stimuli  will  be  exa^Tiined  for  trends  of  change,  noting  that  variations 
caused  by  the  background  noise  should  be  random  in  nature. 

In  a  preliminary  investigation  of  the  EER's  in  Figure  3.2  as 
a  function  of  stimulus  nuirfDer,  little  transient  change  is  discernible. 
The  only  noticeable  diange  is  a  general  decrease  in  amplitude  between 
tlie  first  and  remaining  responses.     The  first  EER  appears  to  contain 
a  low  frequency  coi^onent,  particularly  evident  in  the  latter  portion 
o£  this  response,  which  is  not  as  predominant  in  the  EER's  that 
follav.     In  the  superimposed  EER's  in  Figure  3.3,  this  la^r  frequency 
activity  cabases  the  initial  components  of  tiie  first  response  to  rise 
above  the  other  responses  and  the  later  components  to  descend  below 
the  remaining  responses.     In  addition,  the  EER  to  the  first  stimulus 
contains  a  sustained  positive  conponent  about  300  irsec,  which  does 
not  have  a  counterpart  in  any  of  tiie  following  EER's.     However,  the 
reliability  of  this  component  is  questionable  because  of  the  increase 
in  tlie  ensemble  standard  deviation  in  this  region.     It  therefore 
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appears  that  not  only  does  the  average  response  to  the  first  stimuli 
differ  from  those  that  follow,  but  that  tlie  predominant  component 
(around  the  third  peak  of  Figure  3.3)   of  tiiis  response  is  more 
variant  than  any  other. 

If  the  activity  outside  of  the  primary  response  region  (60  to 
300  msec)   in  Figure  3.3  is  ignored  and  the  means  within  the  correlated 
section  are  removed  from  eadi  response,  the  amplitude  of  the  first 
response  in  this  region  does  not  stand  out.     For  example,   if  the 
magnitude  of  the  vvaveforms  in  the  primary  response  region  are  plotted 
as  a  function  of  stimulus  number,   the  amplitude  of  the  first  response 
is  no  longer  so  predominant.     A  plot  of  both  the  rms  and  peak  to 
peak,  ptp,   amplitude  variations  for  tlie  individual  response  in  the 
interval  is  provided  in  Figure  3.9.     For  this  figure  the  rms  value  is 
about  a  zero  mean  and  the  ptp  value  was  measured  from  the  second  and 
third  primary  deflections  indicated  in  Figure  3.3. 

Since  the  responses  are  veiy  similar  in  this  region,  the  mis  and 
ptp  measures  provide  essentially  the  same  information.     Both  demon- 
strate a  slight  decreasing  trend  in  amplitude  with  the  number  of 
stimulus  presentations.     However,  the  magnitude  of  the  amplitude 
fluctuations  is  not  much  smaller  tlian  the  overall  decrease;   therefore 
the  significance  of  tlie  latter  appears  questionable.     In  a  plot 
(Figure  3.10)   of  the  clianges  in  the  time  delay  or  latency  of  the 
three  peaks  in  the  60  to  600  millisecond  inter\''al  no  significant 
trend  is  discernible.     Therefore,   the  above  results  imply  that  little, 
if  any,   transient  diange  exists  in  the  first  fifty  responses  for 
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responses  averaged  across  a  subject  population.  An  exception  to  this 
observation  is  the  difference  in  the  overall  character  between  the 
first  and  remaining  responses.  Tliis  may  be  evidence  of  some  initiali- 
zation of  the  visual  system  either  at  the  receptor  or  higher  neural 
processing  level.  At  the  receptor  level,  Brazier  (1967)  has  noted 
that  there  is  a  large  difference  in  the  amplitudes  of  the  electro- 
retinogram  attending  the  first  and  second  stimuli  in  a  repetitive 
train.  However,  she  also  notes  that  the  amplitude  of  the  second  ERG 
is  smaller  than  any  that  follow,  indicating  that  the  transient  change 
at  the  receptor  input  exists  beyond  the  first  flash.  The  persistent 
transient  variations  (11  to  20  seconds)  noted  by  Regan  (]956b)  in 
responses  recorded  at  the  scalp  also  inply  a  more  sustained  initiali- 
zation period  than  is  found  in  the  above  EER  results.  Ho^^?ever,  Regan 
utilized  a  continuous  sinusoidal  stimulus  which  restricts  extra- 
polation of  his  results  to  the  present  study.  If  initialization  is 
not  a  factor  then  a  high  degree  of  variability  in  the  first  response 
might  account  for  the  relative  uniqueness  of  the  EER  to  the  first 
flash.  The  ensemble  standard  deviation  following  the  first  stimulus 
partially  supports  this  reasoning  but  it  does  not  explain  all  of  the 
component  differences  between  the  first  and  succeeding  EER's.  There- 
fore, system  initialization  and  inter-subject  response  deviation  may 
both  contribute  to  this  difference. 

3.7  Average  of  Ensemble  Evoked  Responses 

Tne  results  of  the  previous  sections  indicate  that  the  EER  provides 
a  relatively  stable  response  estimate  in  tiie  region  from  60  to  300 
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msec.     If  this  region  is  stipulated  to  contain  tlie  primary  activity 
of  tlie  evoked  response,   then  very  little  response  information  should 
be  lost  by  averaging  EER  estimates  across  stimuli.     Tliis  average 
should  further  reduce  the  average  background  noise  activity  and 
increase  the  quality  of  the  response  estimate.     That  is,  the  variance 
of  the  background  noise  and  the  variance  of  the  EER  estimate  should 
both  decrease  with  the  increase  in  the  total  number  of  averages 
from  fifty  to  five  hundred. 

Figure  3.11  provides  the  average  of  the  ten  EER's  and  the  ten 
control  waveforms  shown  previously.     As  can  be  seen  in  this  figure, 
this  average  seems  to  extract  the  near  sinusoidal  complex  in  the 
primary  region  and  reduce  the  remaining  components .     Tlie  averaged 
control  is  a  general  overall  reduction  of  the  individual  control  wave- 
forms presented  earlier.     The  S/N  ratio  in  the  primary  response 
region  of  the  averaged  EER  is  estimated  from  the  averaged  control  to 
be  about  5.6.     This  is  a  two-fold  improvement  over  the  average  S/N 
ratio  estimated  for  the  EER  to  a  single  response  from  the  individual 
control  waveforms.     If  confidence  intervals  were  provided  for  the 
averaged  EER, each  interval  would  be  less  than  0.5  yv  from  this 
average.     Therefore  the  average  of  the  ten  EER's  provides  a  reliable, 
relatively  noise- free  inter-subject  response  norm. 

3.8     Summary 

The  average  over  an  ensemble  of  subject  EEC's  (EER)  extracts 
a  relatively  consistent  estimate  of  the  response  to  visual  stimuli. 
The  consistency  of  ten  EER's  is  found  to  be  greatest  for  corrponents 
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in  the  60  to  300  :iisec  region  following  stimulus  onset.  The  variations 
of  these  coiiponents  with  stimulus  nujiiber  are  generally  insignificant, 
indicating  little  cliange,  either  habitual  or  othenvise,  in  the  activity 
evoked  by  the  first  fifty  stimuli.  However,  a  general  modification 
in  the  character  of  the  entire  average  waveform  is  noted  following  the 
EER  to  the  first  stimulus.   It  is  suggested  tliat  this  change  is  the 
result  of  initialization  of  the  visual  system  and  a  greater  variability 
of  the  first  response  across  the  subject  population.  Tliis  variability 
of  the  first  response  is  demonstrated  in  confidence  intervals  about 
the  ten  EER  estimates. 

The  amount  of  noise  activity  in  the  more  consistent  region  of  the 
EER's  is  estimated  by  ensemble  control  averages.  From  these  estimates 
the  average  S/N  ratio  in  tlie  EER  is  found  to  be  2.6.  In  order  to 
further  extract  tlie  response  signal  from  the  noise,  the  average  of 
the  ten  EER's  is  calculated.  This  averaged  EER  provides  a  more 
reliable  response  estimate  than  single  EER's  and  improves  the  S/N 
ratio.  The  response  estimate  extracted  by  this  average  is  almost  an 
inverted  sinusoid,  beginning  approximately  100  msec  after  stimulus 
onset. 


CHAPTER  4 
THEORETICAL  ASSUMPTIONS  UNDERLYING  RESPONSE  AVERAGING 

In  order  to  theoretically  justify  the  results  obtained  from  many 
EEG  data  processing  techniques,  it  is  necessary  to  assume  certain 
hypotheses  about  the  statistical  character  of  the  fluctuating  poten- 
tials.  Investigators  who  employ  the  averaging  techniques  to  extract 
time  locked  signals  from  the  incoherent  background  noise  generally 
assume  that:  (1)  tlie  noise  is  a  wide  sense  stationary  random  process; 
(2)  the  signal  and  noise  are  independent;  and  (3)  the  signal  and 
noise  are  additive.  The  latter  assumption  is  probably  the  most 
inportant  since  non- additive  components  cannot  be  separated  by  the 
linear  process  of  averaging.  However,  a  signal  dependent  on  non- 
stationary  noise  can  be  extracted  if  the  degree  of  dependency  and 
non-stationarity  is  minor. 

In  this  chapter  an  attempt  will  be  made  to  gain  positive  or 
negative  evidence  about  these  assumptions  by  examining  the  data 
obtained  from  the  subject  ensemble.  The  results  of  these  examinations 
will  not,  however,  enable  the  experimenter  to  negate  or  to  confirm 
these  assuiiptions  because  of  the  lack  of  knowledge  about  the  poten- 
tials being  tested.  Tlius,  to  examine  one  of  the  above  assumptions, 
the  test  utilized  must  hypothesize  the  nature  of  the  signal  and 
noise  parameters.   In  many  cases  the  test  hypotheses  depend  on  one 
or  two  of  the  above  assimptions ,  conditioning  the  conclusiveness  of 
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the  result.  The  present  study  is  tlierefore  intended  to  provide 
insight  into  the  possible  workability,  rather  than  ultimate  veracity 
of  tliese  assumptions.  Sucli  an  examination  may  then  provide  a  foun- 
dation for  further  studies  into  the  statistical  properties  of  evoked 
signals  in  background  noise. 

4.1  Stationarity  of  the  Background  Noise 

A  random  process  is  said  to  be  stationary  if  the  statistics  of 
the  process  are  time  invariant  (see  Davenport  and  Root,  1958).  For 
example,  a  large  number  of  statistically  independent  sample  functions 
which  are  simultaneously  recorded  constitute  a  wide  sense  stationary 
process  if  the  mean  and  autocorrelation  function  at  one  instant  of 
time  equals  those  at  all  other  times.  In  the  present  study  the 
sajnple  functions  are  the  EEG  records  of  the  fifty  subjects.  The 
time  varying  function  of  particular  interest  is  the  ongoing  activity 
under  stimulus  conditions.  Havever,  these  potentials  cannot  be 
separated  from  the  evoked  response.  Therefore,  this  investigation 
must  rely  on  the  resting  potentials  under  non-stimulus  conditions  to 
estimate  the  background  activity.  It  then  becomes  necessary  to 
work  under  tlie  constraint  of  the  second  assumption  (i.e.,  that  the 
signal  and  noise  are  independent)  to  examine  the  first. 

Since  an  average  is  an  estimate  of  the  first  moment  of  the 
ensemble,  it  is  the  stationarity  of  background  noise  mean  which  is 
the  most  critical  to  the  averaging  teclinique.  If  tiie  noise  mean  is 
highly  non- stationary  (e.g.,  the  averaged  potential  fluctuations 
exceed  those  of  tJie  response  average)  then  the  averaging  process 
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cannot  extract  a  detectable  signal  from  the  background  noise. 
Utilizing  the  resting  EEG  potentials  to  approximate  the  background 
activity,  tlie  noise  mean  estimate  (ensemble  control  average) ,  seen  in 
Figure  3.4,  is  not  time  invariant.  However,  are  these  fluctuations 
tiie  result  of  an  underlying  non -stationary  process  or  caused  by  the 
variations  inherent  in  the  estimation  process?  In  order  to  answer 
this  question  it  is  necessary  to  determine  the  inherent  variance  of 
the  mean  estimate  of  N  sample  functions.  If  the  N  sample  functions 
are  independent,  which  is  plausible  since  individual  subject  EEC's 
should  not  depend  on  one  another,  then  this  variance  of  the  estimate 
is  (see  Bendat  and  Piersol,  1966) 

VAR[m^(t)]=  ^j^  (4.1) 

where  ai(t)  is  the  variance  of  the  process  which  is  unknown  and 
must  be  estimated.  This  latter  estimate  also  has  variance  associated 
with  it  which  again  must  be  estimated,  and  so  forth.  In  order  to 
prevent  endless  estimation  of  these  statistics,  the  estimate  of  the 
variance,  ai(t) ,  will  be  assumed  equal  to  the  actual  variance, 
a|(t),  i.e., 

a|(t)  =  a|(t)  (4.2) 

Now  utilizing  the  Tchebydieff  inequality,  the  range  of  the  mean 
estimate  can  be  predicted  witliout  knowing  its  exact  probability 
distribution.     T!us  inequality  states  that  the  probability  that  the 
difference  between  the  true  mean  and  estimate  mean  will  be  greater 
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than  C  standard  deviations  of  the  estimate  will  be  less  than  1/C^, 
i.e. , 

PROB   [|m^(t)-m^(t)|>  C{sd[ni^Ct)]}]   >  -  (4.3) 

In  Chapter  3  this  inequality  was  used  to  provide  89  per  cent  confi- 
dence limits   (C=3)   about  the  EER  and  \'ER  estimate,  noting  that 

a.(t) 
sd[m  (t)]   =  -^ (4.4) 

Similar  confidence  limits  can  be  obtained  by  hypothesizing  a 
mean  and  using  Or(,t)   to  obtain  confidence  limits.  For  example,  if 
the  noise  mean,  m^(t) ,  is  h>^othesized  to  be  stationary  and  equal  to 
zero,  then  Equation  4.3  becomes 

PROB[|m.(t)|>  C{a^(t)/v^}]  <  i  (4.5) 

t         -  ±  -  ^2 

Again  setting  C=3,  the  89  per  cent  confidence  intervals  about  the 
hypothesized  stationary  mean  can  be  established.  In  addition,  from 
Equation  4.5  another  value  of  C  can  be  calculated  from  the  ensemble 
value  of  mj-(t) ,  the  ensemble  control  average,  and  ar(t) ,  the  control 
standard  deviation,  i.e., 

m.(t) 

C  <  v^  — (4.6) 

a^(t) 

The  maximum  value  of  C  found  in  Equation  4.6  can  then  be  compared 
to  the  tlireshold  value  of  3  --  rejecting  the  stationary  hypothesis 
for  C>3  and  accepting  the  hypothesis  for  C<3.  This  is  equivalent 
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TABLE  4.1 
Maximum  Number  of  Standard  Deviations  of 
Ensemble  Control  Average  About  a  Zero  >fean 

Control     1     2     3    4     5    10    20    30    40    50 
C^^^      2.56  1.62  2.53  2.43  2.24  3.00  2.37  3.07  3.69  2.37 

IllaX 

to  testing  the  h)T)othesis  at  a  significance  level  of  greater  than 
5.5  per  cent  since  the  probability  of  acceptance  of  this  two- tailed 
test  is  greater  than  1-1/C^,  if  the  hypothesis  is  true.  If  the  under- 
lying distribution  is  Gaussian,  this  significance  becomes  greater 
than  1  per  cent.  The  maximum  value  of  C  for  each  of  the  control 
trials,  whose  waveform  is  proAdded  in  Figure  3.4  is  given  in 
Table  4.1.  From  this  table  the  stationary  mean  hypothesis  is 
accepted  seven  times,  rejected  twice,  with  one  indecision.  These 
results  would  therefore  iiiply  that  the  ensemble  mean  of  the  resting 
EEC,  and  perhaps  the  background  activity,  is  nearly  time  invariant 
with  an  occasional  non-stationar)'-  trend. 

Another  estimate  of  the  degree  of  fluctuation  of  the  first 
moment  of  the  background  activity  is  the  plus -minus  average  suggested 
by  Schimmel  (1967) .  This  estimate  is  similar  to  the  HER  except  that 
alternate  saiTiple  fonctions  are  subtracted  in  the  averaging  process 
rather  than  added.  This  estijnate  can  be  expressed  mathematically  by 

-j   odd  even 

±nyt)   =  ^1  I  f.(t)   -     Z     f.(t)]  (4.7) 

1  3       ^ 

If  the  evoked  response  is  assumed  to  be  nearly  the  same  for  each 
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subject,  i.e.,  near  deterministic,  then  the  resultant  average  should 
contain  mainly  background  activity.  The  reliabilit)'  of  this  estimate 
is,  however,  contingent  upon  the  third  assumption  (i.e.,  that  the 
signal  and  noise  are  additive)  in  addition  to  the  second. 

With  this  restriction  in  mind,  the  preceding  test  for  station- 
arity  can  also  be  applied  to  the  plus -minus  average.  The  standard 
deviation  of  this  estimate  about  the  true  +_  average  is  the  same  as 
the  standard  deviation  of  mr(t)  provided  in  Equation  4.4  since  ;Kj|(t) 
is  the  same  as  the  variance,  ar-(t) ,  of  the  process  (see  Schiramel) 
under  stimulus  conditions.  By  assuming  the  response  signal  to  be 
near  deterministic,  this  variance  is  mainly  the  result  of  the  back- 
ground activity.  If  we  again  assume  that  the  estimate  of  the  variance 
is  very  near  the  actual  variance,  the  values  for  C    can  be  calculated. 

The  values  for  C    provided  in  Table  4.2  again  indicate  that  the 

max  ^ 

first  moment  of  the  background  activity  is  stationary,  i.e.,  the 
stationary  hypothesis  is  accepted  for  every  +_  average. 

If  the  signal  is  assumed  to  be  non- deterministic,  it  contributes 

to  the  estimated  variance  of  the  process  and  the  below  C    values 

TABLE  4.2 
Maximum  Number  of  Standard  Deviations 
of  +_  Estimate  about  Zero  Mean 

Response     1     2     3     4     5    10    20    30    40    50 

C        1.93  2.34  2.69  2.00  2.21  2.14  2.30  2.49  2.27  2.82 
max 
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increase.     That  is,  tiie  noise  variance  is  now  the  estimated  ensemble 

variance,  a?.(t) ,  minus  the  signal  variance.     In  addition,   the  ±  average 

will  contain  some  signal  activity.     In  order  to  accept  the  stationary 

mean  hypothesis  for  the  ten  ±  estimates  tlie  standard  deviation  of  the 

signal  cannot  be  greater  than  10  per  cent  of  the  total,  aj-(t) .     Such  a 

limitation  is  reasonable  since  present  magnitude  estimates  of  the 

evoked  response  indicate  that  the  background  EEG  activity  has  an 

amplitude  about  five  to  ten  times  as  great  as  the  response   (Kitasato 

and  Hatsuda,   1965;  Perry  and  Childers,  1969).     Therefore,  with  the 

incoherent  nature  of  the  noise  and  time  locked  character  of  the 

signal,  the  ensemble  variability  of  the  foraier  should  be  at  least 

ten  times  greater  tlian  the  variations  of  the  latter. 

In  order  to  check  tlie  validity  of  the  previous  test,  it  is 

interesting  to  examine  the  stationarity  of  the  response  average,  i.e., 

EER.     Again,   assuming  that  the  signal  and  noise  are  additive  and 

independent,  and  that  the  signal  standard  deviation  is  at  most  10  per 

cent  of  the  total  ensemble  variance,  then  the  values  of  C         for  the 

'  max 

ten  EER's  in  Figure  3.2  are  given  in  Table  4.3. 

TABLE  4.3 

Maximum  Number  of  Standard  Deviations 

(times  0.1)  of  EER  Estimate  About  Zero  Mean 

Response     1     2     3    4     5    10    20    30    40    50 

C         42.5  36.1  28.1  30.1  22.6  38.3  33.0  35.8  28.6  19.0 
max 
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These  results  clearly  demonstrate  that  the  mean  of  the  signal  with 
noise  appears  to  be  non- stationary,  which  is  imperative  for  the  signal 
to  be  extracted  by  averaging.  Even  if  the  signal  variance  accounts  for 
the  entire  ensemble  variance,  the  C^^  values  calculated  (i.e.,  1/10 
of  the  above  results)  would  indicate  the  EER  to  be  non- stationary 
six  out  of  ten  times.  Therefore,  the  preceding  stationary  mean  test 
can  detect  non- stationary  trends  under  the  worst  case  conditions, 
which  indicates  that  it  is  a  valid  test. 

In  order  for  the  background  activity  to  be  wide  sense  stationary 
the  autocorrelation  function  must  also  be  time  invariant.  The  auto- 
correlation function,  R(t,  t-x) ,  is  a  function  of  both  time,  t, 
and  the  time  shift  of  the  correlation,  x.  For  R(t,  t-x)  to  be  time 
invariant  the  autocorrelation  function  at  time  t^  must  be  the  same  at 
all  other  time  \^   for  a  given  x,  i.e., 

RfCtp  t^-x)  =  Rf(t2,  t2-T)  C4.8) 

The  estimate  of  this  function  is  obtained  utilizing  Equation  4.9. 

Let,  t-x)  =  \    Z  fi(t)fi(t-x)  (4.9) 

i=l 

Using  Equation  4.9,  a  three-dimensional  plot  of  the  ensemble  auto- 
correlation function  estimate  versus  t  and  t-x  for  the  resting  EEG 
ensemble  (Control  1)  is  provided  in  Figure  4.1.  Since  the  previous 
test  indicated  that  the  mean  is  stationary,  this  average  is  sub- 
tracted from  the  autocorrelation  to  eliminate  any  sampling 
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fluctuations  of  tiiis  first  moment.  Therefore  tlie  autocorrelation 
function  in  Figure  4.1  is  more  correctly  an  autocovariance  function 
which  has  the  form 

C.(t,  t-T)  =  ^  S  [f.(t)  -  m.(t)][f.(t-T)  -  in.(t-T)]       (4.10) 
•'•  i 

The  magnitude  of  this  function,  like  the  autocorrelation  function,  is 
a  maximum  for  t=0,  which  occurs  along  the  main  diagonal  of  Figure  4.1. 
It  also  should  be  noted  that  the  covariance  function  is  symmetrical 
about  T=0;  therefore  the  super  and  subdiagonal  values  are  equivalent. 
In  order  to  examine  the  stationarity  of  this  function  one  must 
examine  the  fluctuations  of  C^(t,  t-x)  for  a  fixed  x  and  varying  t, 
which  is  analogous  to  viewing  the  variation  C  along  the  main  and  off 
diagonal  terms  in  Figure  4.1.  Fluctuations  of  C^(t,  t-x)  with  time 
are  noticeable  in  this  figure;  however,  their  significance  is 
questionable.  From  the  actual  values  of  C^(t,  t-x)  provided  in 
Table  B.21  of  Appendix  B,  one  does  notice  that  the  magnitude  of  these 
fluctuations  is  generally  independent  of  the  diagonal  used,  i.e., 
the  time  shift,  x.  Therefore  fluctuations  of  the  main  diagonal 
for  x=0,  which  is  equal  to  the  ensemble  variance,  should  provide  as 
good  a  measure  as  any  other  C^(t,  t-x^)  of  the  time  fluctuation  of 
the  correlation  function. 

The  ensemble  standard  deviations  of  the  ten  control  trials  are 
superimposed  in  Figure  4.2.  Iliese  deviations  are  overlapped  to  pro- 
vide a  simultaneous  estimate  of  the  degree  of  fluctuation  and  to 
provide  a  comparison  with  the  stajidard  deviations  presented  in  the 
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same  manner  in  Figure  3.6.  Tlie  reader  interested  in  the  individual 
ensemble  standard  deviations  of  eadi  control  is  referred  to  Figure  B.5. 
The  maximum,  minimum,  and  average  values  of  the  waveforms  in 
Figure  4.2  are  4.5  yv,  1.8  yv,  and  2.8  yv,  respectively.  If  the 
standard  deviation  is  hypothesized  to  be  stationary  and  equal  to  the 
average  of  these  waveforms,  tlien  the  fluctuations  of  the  estimated 
values  exceed  a  range  of  61  to  160  per  cent  of  this  average. 
This  variation  is  probably  greater  than  one  would  expect  from  the 
estimation  error;  however,  without  knowledge  of  the  higher  moments 
one  cannot  further  quantify  this  statement. 

These  higher  moments  are  essentially  descriptors  of  the  finer 
detail  of  the  underlying  probability  structure  of  the  process. 
Therefore  the  profile  variations  of  the  ensemble  probability  density 
function,  pdf ,  should  provide  a  rough  evaluation  of  the  degree  o£ 


1  -^ 


FIGURE  4.2 

Superimposed  Ensemble  Standard  Deviations 

Each  standard  deviation  attends  one  of  the  ten 
occluded  stinuali  which  occur  at  the  extreme 
left  of  the  figure. 
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time  dependence  of  these  moments.  Previous  investigations  o£  the 
amplitude  probability  density  function  of  individual  EEG  records 
have  provided  conflicting  results.  Elul  (1968)  and  Saunders  (1963) 
note  that  EEG  amplitude  histograms  appear  to  approximate  the  normal 
(Gaussian)  distribution.  However,  the  results  of  a  study  by  Campbell, 
et  al.  (1967)  generally  reject  the  hypothesis  that  the  EEG  ajiplitude 
is  noiTTially  distributed.  A  comparison  of  tlie  ensemble  histograms  to 
the  normal  distribution  may  provide  additional  information  about  the 
underlying  EEG  pdf . 

Histogram  estimates  of  the  ensemble  pdf  at  every  10  millisecond 
interval  after  the  first  occluded  stimulus  are  sho\vn  in  Figure  4.3. 
Eadi  of  the  seven  ravs  in  this  figure  approximately  correspond  to 
the  seven  135  millisecond  intervals  used  previously  to  segment 
the  ensemble  control  averages  (Figure  3.4).  As  can  be  seen  in 
Figure  4.3,  there  appears  to  be  little  consistency  between  the 
ensemble  histograms  at  different  time  instances.  Tae   significance 
of  these  pdf  deviations  with  time  is  conditioned  by  the  variance  of 
tiie  estimation  process  (i.e.,  the  histogram  approximation).  However, 
tlie  wide  contrast  in  contour  between  individual  pdf  estimates  appears 
to  exceed  the  limitation  of  sample  variability.  For  example,  the 
histograms  in  Figure  4.3  vary  from  unimodal  (t=60  msec)  to  bimodal 
(t=90  msec) ,  as>nnmetric  (t=710  msec)  to  near  symmetric  (t^lOO  msec) , 
and  leptokurtic  (t=490  msec)  to  platykurtic  (t=330  msec) . 

In  addition  to  this  variability,  few  of  the  histograms  resemble 
the  normal  distribution.  An  attempt  to  quantify  the  degree  of 
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noraiality  utilizing  the  Kolmogorov-Smimov  test  proved  unsuccessful. 
Although  the  Kolmogorov-Smimov  test  is  generally  more  efficient  than 
the  chi-square  goodness  of  fit  test  (Miller  and  Freund,  1965)  for  the 
number  of  available  samples  (50) ,  the  necessity  of  estimating  the 
mean  and  variance  reduced  the  sensitivity  of  the  test  (see  Lilliefors, 
1967).  That  is,  there  was  a  very  high  probability  of  accepting 
the  hypothesis  of  normality  when  it  was  false. 

4 . 2  Dependence  of  Signal  and  Noise 

In  the  previous  section  we  noted  that  if  the  signal  and  noise 
are  additive  and  the  signal  near  deterministic,  then  the  plus-minus 
average  should  provide  an  average  estimate  of  background  activity. 
If  the  same  averaging  process  is  performed  on  the  resting  EEG, 
then  a  similar  estimate  of  the  noise  under  non-stimulus  conditions 
is  obtained.  A  conparison  of  the  general  character  of  these 
processed  potentials  should  provide  an  estimate  of  variations  in 
the  average  noise  under  stimulus  and  non-stimulus  conditions.  How- 
ever, how  does  one  make  such  a  comparison?  The  waveshape  of  these 
potentials  for  ten  signal  and  ten  control  averages  (see  Figures  B.3 
and  B.4)  has  little  meaning,  since  they  were  recorded  at  different 
times.  But  tlie  degree  of  fluctuation  in  tlie  rms  measure  of  the  ± 
average  should  provide  an  estimate  of  average  energy  under  both 
conditions.  Table  4.4  provides  the  rms  values  of  each  plus  minus 
control  and  response  waveform.  Tlie  range  of  rms  values  is  nearly 
equal  under  both  conditions.  Thus,  tlie  time  fluctuating  averages 
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TABLE  4.4 
R>5  Values  in  yv  of  ±  Averages 
Under  Stimulus  and  Non- Stimulus  Conditions 

Waveform    1     2     3     4     5    10    20    30    40    50 
±  Response  0.46  0.37  0.42  0.28  0.39  0.38  0.32  0.36  0.27  0.46 
±  Control  0.26  0.36  0.40  0.38  0.32  0.28  0.49  0.58  0.45  0.67 

of  the  background  activity  and  resting  activity  are  of  the  same 
magnitude.  Since  the  signal  potentials  are  present  in  one  case  and 
not  the  other,  the  general  implication  is  tliat  the  magnitude  of  the 
average  noise  is  independent  of  the  signal. 

Another  property  of  the  background  and  resting  EEC  which  can 
be  investigated  under  the  assumptions  of  a  near  deterministic  signal 
added  to  noise  (s+n) ,  is  the  ensemble  standard  deviation.  If  the 
signal  and  noise  are  independent,  the  standard  deviation  of  the 
resting  EEC  ensemble  (Figure  4.2)  should  have  a  character  similar 
to  the  standard  deviation  of  the  stimulus  provoked  activity 
(Figure  3.6).  That  is,  the  ensemble  standard  deviation  of  signal 
plus  noise ,  a  ^  (t) ,  is  given  by 


o^^^it)  =   [a|(t)  +  o^(t)  +  2a^^(t)]^/2        ^^-^^^ 

If  the  signal  is  near  detenninistic  and  independent  from  the  noise, 
then  the  signal  variance,  a^(t),  should  be  much  less  than  the  noise 
variance,  a^(t),  and  the  covariance  between  the  signal  and  noise, 
o^    (t) ,  equals  zero,  which  reduces  Equation  4.11  to  the  follaving: 
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a   (t)  -  a  (t)  (4-12) 


Thus,  the  standard  deviation  of  the  signal  plus  noise  ensemble 
is  primarily  the  result  of  the  background  noise  and  should  be  similar 
to  the  standard  deviation  of  the  resting  EEG,  if  the  stimulus  does  not 
affect  the  noise  [i.e.,  signal  and  noise  are  independent).  As  noted 
previously,  the  time  fluctuations  of  the  standard  deviations  should 
not  be  tlie  same  for  both  cases,  since  these  deviations  were  recorded 
at  different  times.  If  we  compare  the  standard  deviations,  a^^^Ct) 
and  a  (t) ,  provided  in  Figures  3.6  and  4.2  respectively,  the  overall 
average  and  peak  to  peak  values  are  veay  similar,  i.e., 

ave[ag^^(t)]  -  2.82  yv       ave[a^(t)]  =  2.86  yv 
ptpLa^^^Ct)]  =  3.1  yv        ptp[a^(t)]  -  2.7  yv 
From  these  measurements  one  would  predict  that  the  variance  of  the 
noise  is  independent  of  the  signal.  However,  in  Figure  3.6  a  slight 
increase  in  the  average  standard  deviation,  5^^^{t) ,   was  noted  around 
the  primary  response  region.  Since  similar  trends  do  not  appear  in 
the  standard  deviations,  a^(t) ,  in  Figure  4.2,  this  increase  could 
imply  some  signal  and  noise  dependence. 

If  the  near  deteministic  signal  and  noise  are  dependent  and  the 
average  value  of  the  noise  is  small,  then  the  ensemble  standard 
deviation  of  signal  plus  noise  becomes 

a       (t)  =  (oHt)  +  E[s(t)n(t)]}  (4.13) 

s+n       n 

where  E[s(t)nCt)]  is  the  mathematical  expectation  or  stochastic 
average  of  the  product  s(t)  n(t) . 


59 


The  ensemble  standard  deviation  for  noise  dependent  on  a  deter- 
ministic signal  should  then  increase  when  the  ainplitude  of  the  signal 
increases.  This  is  the  case  in  Figure  3.6,  particularly  for  the 
standard  deviation  of  the  first  response  (see  Figure  3.1).  However, 
if  the  signal  is  relatively  non-deterministic,  then  the  standard 
deviation  of  independent  signal  and  noise  would  also  tend  to  increase 
in  the  region  of  greatest  signal  variability.  Since  it  is  possible 
that  the  standard  deviation  of  a  random  signal  could  be  proportional 
to  the  magnitude  of  the  signal,  the  existence  of  a  random  signal  would 
also  explain  the  slight  increase  in  Figure  3.6.  In  either  case  this 
correspondence  between  signal  magnitude  and  the  ensemble  standard 
deviation  is  minor,  except  for  the  first  response;  tliis  would  imply 
that  the  standard  deviation  of  the  noise  is  relatively  independent 
of  the  signal. 

A  comparison  of  auto co variance  functions  under  stimulus  and 
non- stimulus  conditions  did  not  reveal  any  additional  inform.ation 
about  the  dependence  of  the  signal  and  noise.  The  variations  of  this 
function  appeared  to  depend  primarily  on  the  variance,  a£.(t)  ,  as  noted 
previously  (see  Figures  4.1,  B.6  and  B.7).  The  histogram  estimates 
of  ajnplitude  probability  density  functions  under  stimulus  and  non- 
stimulus  conditions  were  too  variable  to  attempt  a  comparison  (see 
Figures  4 . 3  and  B . 8) . 

4.3  Additivity  of  Signal  and  Noise 

In  the  previous  chapter  the  maximum  (9.68),  minimum  (1.03)  and 
average  (2.6)  rms  signal  to  noise,  S/N,  ratios  wei^e  computed.  Tliese 
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ratios  were  calculated  from  the  ratio  of  the  mis  value  of  the  EER's 

in  the  second  intei'val  to  the  rms  value  of  the  ensemble  control 

average  in  any  interval.      In  addition,  the  correlation  coefficients 

were  found  bet^veen  EER  wav^eforms.     Since  the  S/N  ratio  estimates 

tlie  ajTiount  of  noise  remaining  in  the  EER  waveforms ,   a  decrease  in 

the  S/N  ratio  should  be  reflected  by  a  decrease  in  the  correlation 

coefficients  between  waveforms.     Tliat  is,  if  tlie  signal  and  noise 

are  additive,  then  the  correlation  between  EER's  should  depend  on 

the  relative  amount  of  uncorrelated  noise  and  correlated  signal 

contained  in  these  averages. 

In  Chapter  3  the  correlation  coefficient  between  EER  waveforms 

X  and  y  was  calculated  from  the  time  standard  deviations,  ¥     and  T  , 
/  '     X  y' 

and  the  time  covariance  T     ,   i.e., 

P=^  (4.14) 

X  y 

The  S/N  ratios  were  found  using  the  equation 

R     ¥ 

S/N  =  -^=J-  (4.15) 

rms  "H 
n 

where  H'  is  an  estimate  of  tlie  average  signal  time  standard  deviation 
and  H*  is  an  estimate  of  the  average  noise  deviation  under  stimulus 
conditions  (assuming  the  signal  and  noise  are  independent) .  Both  of 
these  functions  are  dependent  on  the  time  variance,  T,  and  a  mathe- 
matical relationship  bef.veen  the  two  can  be  established  with  certain 
constraints. 
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If  the  EER  is  composed  o£  signal  plus  noise,   i.e., 

X  =  Si  +  ni  (4.16) 

y  =  32  +  na  (4.17) 

then  the  time  covariance  is 


^      =  ^  -  Y         -  1}/         +  1}/  (4.18) 

xy        S1S2         Sin2       S2ni       nin2 

If  the  signal  and  noise  are  statistically  independent,  which  insures 

linear  independence,  then  the  tv^'o  covariances  between  the  signal  and 

noise  in  Equation  4.18  are  zero  and  Equation  4,14  becomes 

^        +4/ 
S1S2       nin2 
p  =  ^/2  ^4.19) 

[(^,  +\  )(y^  +\  )r 

Si     ni       S2     n2 

If  the  noise  is  stationaiy  then  the  time  standard  deviation,  ^  ,  of 
its  estimated  average  over  equal  periods  of  time  should  be  the  same, 
providing  the  period  is  sufficiently  long.     If  we  assume  that  the 
period  between  stimuli  is  sufficiently  long  and  that  the  ensemble 
average  signal  is  near  deterministic  for  each  stimulus,  then 


and 


'^Si   ""  '^sz       '^s  (4.20) 


'i'       -  Y       =  "V  (4.21) 

na         n2         n 


Tlie  correlation  coefficient  betiveen  EER  waA^efomis  can  tlien  be  written 
as  in  Equation  4.22. 
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S1S2       nin2 
p  =  C4.22) 

s        n 

A\'hich  reduces  to 

p   (S/N)^  +  p„ 
p  =     s  ^  (4.23) 

1  +   (S/N)^' 
where  p     and  p     are  the  correlation  coefficients  bet^s'een  the  signal 
and  noise  and  S/N  is  defined  by  Equation  4.15  with  the  "  deleted. 
Since  we  assumed  previously  tiiat  the  average  signal  was  near  deter- 
ministic,  p     is  almost  unity.     The  noise  should  be  incoherent  with 
the  stimulus;  tlius  p     is  near  zero  and  Equation  4.23  reduces  to 


(S/N) ^ 

p  =  (4.24) 

1  +   (S/N) 2 

Therefore,  if  a  pair  of  identical  signals  are  added  to  stationary 

noise,  which  is  independent  from  the  signal,  then  the  correlation 

coefficient  between  the  resultant  wavefoims  can  be  calculated  from 

the  S/N  ratio. 

Using  the  three  S/N  ratios  given  on  page  59,   the  following 
values  of  p  are  found  from  Equation  4.24: 

p         =  0.989 

p   .      =  0.515 

p         =0.871 

Therefore,  subject  to  the  above  stipulations,  one  would  predict  the 
above  values  for  correlation  coefficients  between  the  EER  responses 
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in  the  second  interv-al.     A  comparison  of  tiiese  predicted  correlations 

with  those  calculated  in  the  previous  chapter,   i.e., 

p         =  0.99 

p   -     =  0.53 
min 

p         =0.88 

discloses  that  the  former  and  latter  are  nearly  equivalent. 

The  implications  of  this  strong  similarity  between  the  predicted 

and  measured  coefficients  are  two-fold.  First,  the  R    to  C 

'  rms  rms 

signal  to  noise  ratio  provides  a  good  estimate  of  the  background 

activity  remaining  in  the  EER.     Second,   the  assumption  of  an  additive 

signal  to  noise  is  furtlier  supported.     One  argument  which  might  be 

presented  against  the  validity  of  the  predicted  coefficients  is  that 

the  assunptions  of  p  =0  and  p  =1.0  are  not  realistic.     However,  the 
^  n  s  ' 

average  correlation  coefficients  calculated  across  segments  of  the 
control  waveforms  were  very  near  zero.     This  would  partially  reject 
the  argument,  particularly  in  the  case  of  the  average  predicted 
correlation  coefficient. 

Aiiother  method  to  investigate  the  additivity  of  the  signal  and 
noise  is  to  observe  the  growth  of  the  signal  to  noise  ratio  with  the 
number  of  averages.     For  example,   if  a  near  deterministic  signal  and 
random  background  noise  are  additive  and  independent,  then  the 
waveform  of  the  average  signal  will  remain  essentially  constant  and 
the  variance  of  the  average  noise  will  decrease  directly  with  the 
number  of  simiJTiati ons  in  the  average.     If  we  assume  tiiat  the  noise  is 
stationary,   tlien  the  true  noise  mean  is  a  constant,  which  in  the 
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present  study  is  preset  to  zero.     ITie  ms  value  of  the  average  noise 
will  approach  this  zero  value  as  the  nuniber  of  smiraations,  N,  becoines 
very  large   (noting  that  tiie  variance  of  this  average  about  tiie  zero 
niean  decreases  by  a  factor  of  i) .     Hie  ™s  value  of  the  average  signal 
should  be  relatively  invariant  witii  the  number  N.     Therefore  the 
sional  to  noise  ratio,  S/N,  should  be  expected  to  increase  as  the 
number  of  sulfations  in  the  average  increases. 

The  minimum  rate  of  the  S/N  ratio  increase  can  be  specified  in 
tenns  of  the  true  ensen^le  variance,  a|(t) ,  and  a  probability  deter- 
mined from  the  Tchebycheff  inequality   (see  Equation  4.3).     Ihus ,  we 
should  expect  the  estimated  S/N  ratio  to  reflect  a  random  increasing 
trend  with  the  nun^er  of  summations  in  the  average.     However,  how 
can  this  increase  be  examned?     The  technique  chosen  to  observe  this 
growth  was  to  calculate  the  rms  value   (in  the  60  to  300  mllisecond 
region)   of  the  EER  at  every  even  nujnber  of  summations  and  to  compare 
this  value  to  the  mis  value  of  the  ±  average  of  the  same  data.     After 
a  number  of  summations ,  the  ±  average  of  the  data  should  provide  an 
amplitude  estimate  for  the  background  noise  and  the  EER  a  signal 
amplitude  estimate.     Tlie  S/N  ratio  can  then  be  formed  by  dividing  the 
EER  by  the  ±  average  at  every  even  number  of  summations,  N,  within 
the  average.     Hiis  procedure  was  repeated  for  each  set  of  data 
following  the  ten  selected  stimuli. 

A  plot  of  the  S/N  ratio  versus  the  number  of  summations  N  for 
the  EER's  to  the  first  five  stimuli  is  provided  in  Figure  4.4.  Althou^ 
the  S/N  ratio  demonstrates  an  overall  slight  increase  with  N  in  each 
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FIGURE  4.4 

S/N  Ratios  Versus  Number  of 
Summations 

S/N  Ratios  for  the  EER  are 
calculated  at  every  even 
number  of  summations ,  N. 
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of  tlie  plots  in  this  figure,  the  random  fluctuations  of  S/N  over- 
shadow this  growth.   In  addition,  it  should  be  noted  that  the  S/N 
ratio  fluctuations  of  the  first  response  are  nearly  twice  those  of  the 
remaining  responses.  In  a  similar  plot  of  the  S/N  ratio  for  EER's  to 
later  stimuli  (see  Figure  B.9),  an  overall  decrease  in  this  ratio  with 
N  can  be  noted.  These  plots  appear  to  discredit  the  signal  plus 
noise  assumption.  However,  this  conclusion  contradicts  the  previous 
results  in  this  chapter  and  questions  the  results  achieved  by  the 
established  technique  of  response  averaging.  Therefore,  the  adequacy 
of  the  above  test  becomes  liighly  suspect. 

The  inadequacy  of  this  test  appears  to  be  the  result  of:   (1)  the 
large  variations  in  amplitude  of  individual  EEC's  which  appear  to 
introduce  abrupt  changes  in  the  noise  estimates;  and  (2)  the  S/N 
ratio  estimate  being  probably  inaccurate  except  for  a  large  number 
of  sums  (i.e.,  EER  is  more  noise  than  signal).  The  first  short- 
coming might  be  eliminated  by  performing  the  same  test  on  a  single 
subject  (i.e.,  using  sequential  as  opposed  to  ensemble  averaging); 
the  latter  shortcoming  could  be  reduced  by  beginning  with  forty  to 
fifty  sums  and  increasing  the  number  of  summations.  Still  another 
technique  would  be  to  utilize  Equation  4.24  and  discuss  the  improve- 
ment in  S/N  by  the  degree  of  correlation  of  the  average  after  N 
summations  with  a  preset  response  template  or  norm.  Tliis  correlation 
teclinique  would  be  more  sensitive  to  the  waveshape  of  the  response, 
which  would  help  to  separate  signal  and  noise  energy  in  the  case 
of  fewer  summations. 
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4.4  Sununary 

Tlie  separation  of  the  foregoing  examinations  into  three  groups 
is  somev.hat  artificial  because  of  the  inter- dependence  of  the 
results  obtained.  Tnerefore,  many  of  the  results  obtained  actually 
provide  tests  for  more  than  a  single  assumption.  Tlie  above  examina- 
tions are  grouped  into  three  areas  to  help  organize  the  test  results 
and  produce  a  more  readable  text. 

TJie  results  of  these  tests  strongly  imply  that  the  mean  of  the 
background  noise  is  time  invariant  and  that  its  amplitude  is  rela- 
tively independent  of  tlae  signal  activity.  The  higher  order  statistics 
(i.e.,  standard  deviation,  autocovariance  and  pdf)  of  the  noise  appear 
to  be  non-stationary.  However,  the  time  fluctuations  of  the  estimates 
of  these  higher  order  statistics  are  dependent  on  the  variance  of  the 
estimation  process,  which  was  not  investigated.  The  average  and 
peak-to-peak  fluctuations  of  the  ensemble  standard  deviations  under 
stimulus  and  non-stimulus  conditions  were  very  similar,  which  further 
indicates  tliat  the  signal  is  independent  of  the  noise.  The  average 
correlation  coefficients  of  a  theoretical  signal  plus  noise  were 
nearly  identical  to  those  calculated  from  the  actual  averaged  res- 
ponses. Tlius,  the  assumption  of  signal  plus  noise  was  partially 
supported.  However,  a  subsequent  test  of  tlie  theoretical  increase 
in  the  S/N  ratio  by  averaging  signal  plus  noise  data  provided  incon- 
clusive results. 


CHAPTER  5 

SUMMARY 

To  date,  most  investigations  of  evoked  potentials  monitored  at 
the  hiunan  scalp  have  utilized  the  average  of  a  series  of  responses 
recorded  from  a  single  subject  to  define  components  of  tlie  evoked 
response.  This  technique  has  provided  a  fairly  consistent  estimate 
of  an  individual's  average  ER.  Havever,  with  this  method  it  is  not 
uncommon  to  obtain  quite  dissimilar  averages  from  different  subjects 
under  identical  experimental  conditions.  This  contrast  in  the  ER 
estimates  across  a  population  has  hindered  the  interpretation  of  the 
effects  of  stimulus  parameters  on  the  ER  and  frustrated  the  attempt 
to  use  this  average  for  clinical  diagnosis.  It  may  be  possible  that 
the  character  of  an  ER  is,  like  a  signature,  unique  for  each  subject. 
If  this  analogy  is  valid,  utilization  of  the  ER  as  a  research  or 
diagnostic  tool  might  be  limited.  However,  if  one  could  establish  a 
particular  feature  of  ER  estimates  common  to  a  considered  normal 
subject  population,  then  tliis  template  or  nonn  might  identify  that 
quality  associated  with  a  population  of  responses  evoked  by  a  specific 
stimulus.  This  would  be  analogous  to  identifying  the  diaracteristics 
attributed  to  a  specific  letter  of  the  alphabet  in  a  set  of  signatures. 
The  influence  of  variations  in  stimulus  parameters  and  physiological 
abnormalities  might  tlien  be  recognizable  in  tlie  degree  of  deviation 
of  a  subject's  ER  estim.ate  about  this  norm. 


68 


69 


The  present  study  investigated  one  method  of  estimating  such 
a  response  norai.  The  estimation  technique  utilized  was  an  average 
across  an  ensemble  o£  subject  EEC's  containing  evoked  activity.  The 
advantages  of  this  method  were  tv^'o-fold.  First,  it  pemiitted  estima- 
tion of  common  inter-subject  response  activity  evoked  by  a  single 
stimulus.  Second,  the   collection  of  subject  EEC's  strongly  resembled 
an  ensemble  of  independent  sample  functions  from  a  stochastic  process. 
Tliis  second  advantage  permitted  an  investigation  of  several  prevalent 
assumptions  about  the  statistical  nature  of  the  evoked  signal  and 
background  noise. 

The  ensemble  average  revealed  a  fairly  consistent,  recurring 
component  complex  following  each  of  a  series  of  photic  stimuli,  i.e., 
ten  stimuli  diosen  between  the  first  and  fiftieth  stimuli  in  a 
periodic  stimulus  train.  This  averaged  complex  or  primary  response 
mainly  consisted  of  a  positive  to  negative  to  positive  potential, 
near  sinusoidal  in  form,  which  began  approximately  60  milliseconds 
after  the  onset  of  each  stimulus.  Significant  correlation  coefficients 
were  found  bet\veen  all  pairings  of  the  ten  ensemble  averages  in  this 
primar)'  region.  Outside  of  this  region  most  of  the  paired  correlations 
between  waveforms  were  insignificant.  The  average  of  these  paired 
correlations  provided  a  reliable  measure  of  the  conformity  across 
segments  of  the  ten  responses.  Control  estimates  of  the  amount  of 
background  noise  contained  in  tlie  primary  response  region  predicted 
tliat  this  complex  could  contain  from  10  to  50  per  cent  noise,  with 
an  expected  amount  of  28  per  cent. 
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Under  the  assumption  that  tlie  evoked  signal  and  background  noise 
are  additive,  this  prediction  was  furtlier  substantiated  by  tlie  range 
of  the  correlation  coefficients  between  paired  response  estimates. 
Tlierefore,  these  results  implied  tliat  the  control  estimate  provided 
a  good  evaluation  of  the  amount  of  distorting  noise  remaining  in 
the  ensemble  evoked  response,  EER.  Since  the  primary  response  was 
detectable  in  each  EER,  these  results  also  imply  that  even  with 
a  signal  to  noise  ratio  of  unity  the  average  response  is  observable. 

The  ensemble  average  tedinique  also  pro\dded  an  estimate  of 
the  changes  in  the  evoked  response  with  the  number  of  stimulus 
repetitions .  That  is ,  each  EER  was  the  average  of  tlie  activity 
attending  a  single  stimulus.  An  examination  of  the  changes  in  the 
response  estimates  following  each  stimuli  disclosed  moderate  com- 
ponent fluctuations  in  ajnplitude  and  latency  of  the  primary  response. 
Since  these  fluctuations  did  not  exhibit  any  specific  trends  or 
transient  variations  with  the  number  of  stimulus  repetitions,  they 
appeared  to  be  the  result  of  background  noise  distortion  rather  than 
response  changes.  Tlie  one  noted  exception,  the  uniqueness  of  the 
EER  to  the  first  flash,  is  hy'pothesized  to  be  the  result  of  inter- 
subject  response  variability  or  some  initialization  of  the  visual 
system.  This  first  response  estimate  appeared  to  contain  a  low 
frequenc)'  component  '.vhich  was  not  evident  in  the  remaining  estimates. 
In  addition,  the  time  fluctuations  of  the  ensemble  standard  deviation 
about  the  ensemble  average  resembled  tlie  waveshape  of  tliis  average 
for  the  first  45n  milliseconds  following  tlie  first  stimulus.  This 
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was  a  characteristic  relatively  unique  to  the  standard  deviation 
attending  the  first  stimulus  and  resulted  in  a  decrease  in  the  con- 
fidence of  the  first  EER.  Another  singular  property  of  the  potentials 
attending  the  first  stimulus  was  demonstrated  in  plots  of  the  EER 
signal  to  noise  ratio  as  a  function  of  the  number  of  subjects  in  the 
average.  The  magnitude  of  the  variations  of  the  S/N  ratio  in  the 
plot  for  the  first  response  was  nearly  twice  that  found  in  the 
remaining  graphs.  Thus,  these  graphs  reflected  the  greater  varia- 
bility of  the  initial  potentials.  Because  of  these  unique  properties 
of  the  first  response  estimate  it  may  be  advisable  to  delete  this  res- 
ponse when  comparing  \^R's  across  subjects. 

The  average  of  the  ten  EER's  waveform  appeared  to  provide  a 
reliable,  relatively  noise-free  response  template.  In  this  template 
the  primary  sinusoidal  complex  was  enhanced  and  the  amplitudes  of 
the  remaining  components  were  reduced.  The  character  of  these  reduced 
components  strongly  resembled  that  of  a  similar  single  control 
template.  Thus,  this  response  template  further  supported  the  hypo- 
thesis that  the  common  inter-subject  response  activity  generally 
is  restricted  to  the  primary  response  region.  This  average  response 
template  also  implied  that  investigations  which  compare  the  sequential 
average  obtained  from  different  subjects  should  possibly  restrict 
this  comparison  to  the  60  to  300  milliseocnd  region  after  stimulus 
onset.  However,  since  the  above  template  would  not  extract  inter- 
subject  signal  components  whose  latency  varies  from  one  subject  to 
the  next,  tliis  restriction  might  be  inappropriate  for  some  studies. 
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Three  common  assumptions  about  the  statistical  diaracter  of  the 
signal  and  noise  necessary  to  tiieoretically  justify  response  averaging 
were  examined  using  the  ensemble  of  subject  EEC's.  The  results  of 
these  examinations  generally  supported  the  assumptions  that:  (1)  the 
noise  mean  is  stationary;  (2)  the  signal  and  noise  are  independent; 
and  (3)  scalp  potentials  consist  of  signal  plus  noise.  The  first 
assumption  was  investigated  by  testing  t^vo  ax^plitude  estimates  of 
the  average  background  activity  against  a  hypothesized  stationary 
mean.  Of  the  tu^enty  tijnes  tested,  the  stationary  hypothesis  was 
rejected  only  twice.  The  conclusiveness  of  these  tests  was  dependent 
on  the  second  assumption,  i.e.,  the  independence  of  signal  and  noise. 
The  combination  of  signal  and  noise  activity  is  demonstrated  to  be 
generally  non- stationary.  If  the  signal  and  noise  are  dependent,  the 
hypothesis  tests  on  tlie  +  average  should  demonstrate  this  dependence 
by  rejecting  the  stationary  hypothesis.  Since  this  was  not  the  case, 
the  results  of  the  hypothesis  tests  appear  to  definitely  support  the 
assumption  of  a  stationary  noise  mean.  An  investigation  of  the 
stationarity  of  the  higher  noise  moments  was  generally  unsuccessful. 
Since  the  variance  of  the  estimate  of  these  statistics  was  unknown, 
the  time  dependence  of  these  functions  could  not  be  evaluated.  H^vever, 
tlie  relative  degree  of  the  time  fluctuations  of  the  variance,  auto- 
covariance,  and  pdf  functions  implied  that  the  higher  moments  of  the 
noise  might  be  non- stationary. 

An  investigation  of  the  second  assumption,  the  independence  of 
signal  and  noise,  was  accomplished  by  comparing  the  ensemble  +  average 
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and  standard  deviation  under  stimulus  and  non-stimulus  conditions. 
The  ±  average  provided  an  estimate  o£  the  noise  ainplitude  with  and 
without  tlie  evoked  activity.  The  range  of  amplitude  fluctuations 
of  this  average  was  nearly  equal  under  both  conditions.  The  average 
and  peak-to-peak  fluctuations  of  tJie  erisemble  standard  deviation 
under  stimulus  and  non-stimulus  conditions  were  very  similar. 
Therefore,  both  the  noise  average  and  the  deviation  about  this 
average  appear  to  be  independent  of  the  signal  activity. 

An  examination  of  the  third  assumption,  additivity  of  signal 
and  noise,  was  generally  hindered  by  the  inability  to  separate  these 
two  components.  Tliis  shortcoming  thwarted  an  attempt  to  empirically 
demonstrate  the  theoretical  improvement  of  the  signal  to  noise  ratio 
by  averaging  the  assumed  signal  plus  noise  potentials.  Hoi\rever, 
the  relative  stability  of  the  average  estimates  and  the  previous 
prediction  of  the  signal-noise  ratio  by  the  correlation  coefficient 
between  ^^faveforms  supported  this  signal  plus  noise  assujiption. 

The  results  of  this  investigation  of  an  ensemble  of  subject 
EEC's  which  contained  evoked  responses  have  provided  an  estimate 
of  the  common  inter-subject  response.  The  estim.ated  response  template 
foimd  in  this  study  was  limited  to  responses  evoked  by  visual  stimuli, 
specifically  a  monocular,  periodic  flash  of  light.  The  applicability 
of  this  template  to  other  studies  of  responses  in  this  modality 
is  at  present  uncertain.  However,  the  metiiod  utilized  in  this 
investigation  has  provided  at  least  one  teclmique  for  obtaining  an 
inter-subject  response  norm  to  a  specific  stimulus.  It  is  impractical 
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to  utilize  this  ensemble  technique  in  most  experimental  or  clinical 
investigations  because  of  the  large  subject  population  required. 
Tlierefore,  the  applicability  of  the  technique  is  restricted  to 
isolated  studies.  The  results  achieved  in  sudi  studies,  however, 
should  be  useful  in  many  otlier  investigations  which  compare  the  intra- 
subject  response  average  across  a  subject  population. 

The  ensemble  of  subject  EEC's  appears  to  provide  a  worthwhile 
tool  by  which  to  evaluate  theoretical  assumptions  about  the  statistics 
of  the  EEC.  However,  tlie  degree  of  sample  variation  noted  in  this 
study  seems  to  indicate  that  a  greater  number  of  sample  functions 
would  provide  a  more  reliable  evaluation.  This  is  particularly  true 
in  the  case  of  the  estimation  of  the  higher  ensemble  moments 
and  the  ensemble  probability  density  function.  In  order  to  obtain  a 
significant  improvement,  a  fair  amount  of  increase  in  the   number  of 
sample  functions  would  be  necessar>^  The  advantages  of  such  an 
increase  would  probably  be  oversjiadowed  by  the  difficulties  presented 
in  acquisition  of  this  larger  ensemble. 


APPENDIX  A 
EXPERIME>JTAL  PROCEDURE 

A. 1  Data  Acquisition 

Tlie  data  used  in  this  experiment  was  collected  from  fifty  subjects 
(32  males,   18  females).     The  right  eye  acuities  of  all  subjects  were 
noniial  or  correctable  to  noimal  as  measured  on  a  Snellen  chart.     Sub- 
jects were  asked  their  hand  preference  and  given  a  far  vision  sighting 
test  for  eye  dominance.     One  pair  of  bipolar  electrodes   (Lexington 
C-100)  were  placed  5  cm  apart  on  the  scalp,   located  vertically  on  the 
occiput  midline  with  the  lower  electrode  over  the  inion.     This  con- 
figuration was  chosen  becuase  it  has  been  noted  to  produce  more 
reliable  VER's  across  subject  populations.     Each  subject  was  seated 
upright  in  an  electrically  shielded  Faraday  cage  and  asked  to  concen- 
trate on  the  light  seen  through  the  right  eye-piece  of  a  modified 
troposcope.     During  a  stimulus   (signal)   trial  the  light  consisted 
of  a  xenon  flash  of  intensity  3.6  log  ft-lamberts  subtending  2°15'   in 
an  invariant  incandescent  surround  of  intensity  1.3  log  ft-lamberts 
subtending  18°24'.     The  flash  was  presented  periodically  (approximately 
0.9  Hz)   in  two  trials  and  aperiodically  (average  rate  of  approximately 
0.4  Hz)   in  a  third.     Each  si.gnal  trial  was  followed  by  a  control 
(noise)   trial  identical  to  the  preceding  stimulus  trial  with  the 
xenon  flash  occluded.     After  each  trial  the  subject  was  permitted 
a  rest  period  of  two  to  three  minutes;  total  time  required  of  each 


75 


76 


subject  was  usually  less  than  forty-five  minutes. 

The  potential  difference  monitored  between  the  scalp  electrodes 
was  fed  simultaneously  through  Wo   Grass  P5111  differential  amplifiers 
(channels  1  and  2)  with  a  gain  of  lOOK  and  bandpass  of  0.3  to  50  Hz. 
The  bandpass  of  the  signal  in  channel  2  was  further  restricted  by  a 
Khron-Hite  330B  electronic  filter  from  0.3  to  15  Hz.  TTiese  ainplified 
signals  were  then  fed  into  a  Sanborn  2000  Bl  magnetic  tape  recorder 
and  a  Mnemotron  400C  computer  of  average  transients  (CAT) .  Also 
recorded  on  magnetic  tape  was  the  output  of  a  Grass  S-8  stimulator 
which  provided  an  external  trigger  for  the  CAT  and  a  Grass  PS- 2 
photos timulator  which  powered  the  xenon  flash.  In  order  to  detect 
changes  in  equipment  gain  from  one  subject  to  the  next  a  calibrate 
signal  from  a  Medistor  C-IA  Calibrator  was  processed  in  a  manner 
similar  to  EEG  data  prior  to  the  experim.ental  run  for  each  subject. 

A. 2  Data  Processing 

In  each  trial  the  CAT  sequentially  summed  each  one  second  epoch 
of  the  bioelectric  potentials  following  each  stimulus  trigger.  In 
this  manner  VER  and  control  waveforms  were  provided  for  each  subject. 
These  wavefoms  were  permanently  registered  on  a  Sanborn  701A  Strip 
Chart  Recorder.  From  these  records  the  'T3est"  signal  and  noise  trials 
(under  periodic  stimulus  conditions)  were  chosen.  Judgment  of  "best" 
was  obtained  strictly  on  an  amplitude  basis,  disregarding  the  waveshape 
of  the  sunrned  response  to  avoid  biased  results .  That  is ,  the  larger 
sequential  averaged  signal  and  smaller  sequential  averaged  noise 
amplitudes  were  judged  to  be  the  best. 
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The  unprocessed  analog  signals  from  tliese  better  trials  were 
then  digitized  and  stored  on  a  memory  disc  of  an  IBM  1800  conputer. 
Since  the  disc  could  only  accomodate  a  total  of  approximately  4.4 
(10  )  data  words,  it  became  necessary  to  limit  the  amount  of  analog 
data  stored.  Therefore,  the  storage  of  analog  potentials  (digitized 
at  a  rate  of  200  samples  per  second)  was  restricted  to  ten  near  one 
second  epochs  of  data  for  each  subject  and  each  channel.  The   ten 
epochs  chosen  followed  the  first  through  fifth,  tenth,  twentieth, 
thirtieth,  fortieth  and  fiftieth  stimuli  in  the  periodic  stimulus 
train.  These  ten  epochs  were  chosen  to  provide  a  measure  of  both 
the  initial  (first  few  seconds)  and  prolonged  (forty  to  fifty  seconds) 
clianges  in  tlie  evoked  activity  with  repetitions  of  the  stimulus.  Each 
epoch  consisted  of  the  first  190  digitized  points  following  each  of 
the  ten  stimulus  triggers;  thios  a  single  epoch  is  950  milliseconds 
in  duration. 

In  order  to  eliminate  any  near  DC  potentials  introduced  by  equip- 
ment base  line  shifts,  the  mean  was  removed  from  each  190  point  epoch. 
Equipment  gain  fluctuations  across  the  subject  population  were  also 
eliminated  by  multiplying  the  digitized  potentials  by  a  constant 
inversely  proportional  to  the  magnitude  of  the  recorded  calibrate 
signal.  Computations  on  this  adjusted  array  of  data  provided  most 
of  the  results  contained  in  tliis  study.  Tlie  programs  utilized  to 
obtain  the  data  were  fairly  simple  in  nature  and  are  therefore  omitted 
from  the  text.  Tne   numerical  results  of  tliese  computations  were 
converted  to  graphical  displays  by  a  Conputer  Industries  135  On- Line 
Incremental  Plotter. 
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FIGURE  B.9 

S/N  Ratios  Versus  Number  of 
Summations 

S/N  Ratios  for  tlie  EER  are 
calcxilated  at  every  even 
number  of  summations,  N. 
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1.00 


INTtRVAL   1 
.30 


'••e 


0.57  1.00 

0.7a  0.19  1.00 

0.29  0.61  -0.01.       1.00 

0.67  0.5)  0.76      0.52       1.00 

10      O.J't  -0.55  0.1(7      O.IU       0.C5       1.00 

20      0.22  0.514  -0.51      0.58  -0.17    -0.51.       1. 00 

50      0.56  0.10  0.11*    -0.1.5  -O.UO    -0.55      0.57 

UO      0.85  0.56  0.60    -0.07       0.21       0.00      0.56      0.79       1.00 

50      0.76  0.28  0.58      0.11       0.18       0.29      0.1.9      0.0?      0.89       1.00 


1.00 


RtSP      1  2  5  It  5  10  20  50  liO  50 

1.00 

0.99  1.00  fNTEBVUl   2 

0.90  0.90  1.00  ",„   •   .88 

0.92  0.95  0.97  1.00 

0.77  0.77  0.9J  0.85       I. 00 

10      U.92  0.95  0.92  0.97      0.85       1.00 

20      0.97  0.96  0.88  0.91      0.69      0.88       1.00 

50      0.87  0.91  0.75  3.87      0.58      0.92      0.85      1.00 

..0      0.01  0.90  0.91*  0.95      0.79      O.SC      U.76      0.75       1.00 

50      0.96  0.98  0.88  0.95      0.70      0.95      0.95      0.95      0.90      1.00 


RESP      1  2  5  k              S            10           20           SO           bO           50 
l.OO 

-0.10  1.00  INTERVAL   3 

0.k5  -0.55  1.00  Oj,,  ■    .12 

0.88  0.25  0.58  1.00 

-0.U9  0.7»  -0.71  -0.25       1.00 

10    -0.15  0.5S  -0.21  -0.05       0.68       1.00 

20      0.12  0.55  -0.1.6  0.59      0.19    -0.55       l.On 

50      0.01  O.IU  -0.*7  0.29       0.78       D.-S       0.50       I.JO 

1.0      0.55  0.60  -0.06  0.79      0.55       0.25       0.58       0.81       l.CO 

50    -0.U8  0.11  -0.U8  -0.5*.       J. 02    -6.09       0.05    -l.kk    -0.62       l.no 


INTERVAl.  < 

"...  •  -OS 


"CSV      1              2  5              U              5            10 

1  1.00 

2  -0. 26       1.00 
5      O.fiU    -0.72  1.00 

I*       0.U9    -0.25  0.73  1.00 

5       0.67      0.50  0.1*0  0.66      1.00 

ID       0.55       0.55  O.5..  0.05      0.79       1.00 

20  -0.92      0.51  -0.68  -0.51    -O.H    -0.17       1.00 

50  -0.50    -0.05  -0.^9  -0.91    -0.61*    -0.5U      O.Ob      1.00 

1.0  -0.29    -0.02  0.18  0.56      0.05      0.05       0.k7    -0.75       1.00 

50      0.05       0.15  0.02  «0.2U      0.52       0.88      0.05   -0.  U    -0.20       1.00 


RESP      I 

1.00 


INTERVAL  5 


0.20       1.00 

0. 75    -0.k5  l.CO 

0.01    -0.52  0.5U  1.00 

0.91.       0.55  0.56  -0.15       1.00 

10      0.88       0.5U  0.56  -0.23       0.95       1.00 

20      0.59      0.76  -0.29  -0.51  O.GO       0.7C       1.00 

50    -0.1.8    -0.52  O.OC  0.25  -0.53    -0.76    -0.87       1.00 

1.0      0.58       0.85  -0.29  -0.1*0       0.65       0.71       0.95    -1.71       l.)0 

50    -0.62      0.51  -0.57  -0.02  -0.53    -0.51.    -0.29       0.59    -0.03       1.00 


RtSr      1              2  5  I.              S            IC           20           50           40 

1  1.00 

2  0.86  1.00  IKTERVAl  6 
5  -0.26  -0.51.  1.00  -,„  '  Of 
1*       0.15      0.22  0.52       1.00 

5       9.1.C       0.1*1  -0.21*    -0.47       1.00 

10  -0.82    -0.55  0.52      0.55    -0.53       1.90 

20  -0.04    -0.07  -0.02    -0.65       0.37    -J. 54       l.OO 

50  -0.42    -0.42  0.64       0.77    -0.81       0.70    -0.67      1.00 

40  -0.92    -0.66  0.10    -J.l£    -U.22      0.!5       0.15      0.26      1.00 

50  -0.50    -0.54  0.04       0.55    -0.87       0.62    -0.75      0.75      0.52 


RESP      1             2             3  4             5           10          20           50          40          50 

1.00 

0.57      1.00  INTERVAL  7 

0.07      0.08       l.OO  "jvt  ■  ""' 

-0.15    -0.55       0.40  1.00 

0.40      0.49      0. 55  -0. 22       1.00 

10    -0.15    -0.56       0.41  0.98    -0.29       1.00 

20      0.16      0.17    -0.84  -0.17    -0.26    -0.19       1.00 

30    -0.83    -0.67       0.52  0.42    -0.54       0.44    -0.51       1.00 

IlO    -0.78    -0.55       0.53  0.14       0.08      0.14    -0.65       0.87       1.00 

50    -0.18    -0.25    -0.40  -0.57    -0.44    -0.52    -0.05       0.06      0.11       1.00 


TABLES  B.l  through  B.7 

Response  Correlation  Coefficients 

Tlie  above  tables  provide  the  cor- 
relation coefficients  bet\\reen  the 
ten  response  waveforais  in  Figure 
3.2  for  each  135  msec  interval. 


89 


INTERVAL  2 
0         ■    .80 


fiiP  1  2  }  >i  5  10  20  30  kO  50 

1  l.OU 

i  0.97  1.00 

3  0.86  0.8^  1.00 

k  0.<8  0.90  0.90  1.00 

i  0.73  U.7i  0.87  O.si       l.OU 

lu  0.8d  0.93  0.86  0.87       0.83       1.00 

io  0.90  u.ee  u.77  o.eu     0.39     0.79     1.00 

30  0.8k  0.83  0.71  U.7B       0.63       0.8^       0.78       l.OU 

kO  0.89  0.87  U.9U  0.91       0.71       0.77       0.87       0.65       1.00 

50  0.88  0.8b  0.83  0.89       0.60       0.E5      0.83       0.90       0.78       1,00 


IHTlRV/U.  3 


RESP      1              i  3  b  5            10            20            SO            kO           SO 

1  1.00 

2  -0.35       1.00 

3  U.S7  -0.60  1.00 
k  0.63  -0.36  0.39       1 .  OO 

5  -O.kO      0.1k  -0.30    -0.k5       1,00 

lU  -0.01      0.k3  -0.23    -0.15      0.09       1.00 

20  -0.12       0.1k  -0.25    -0.06      0.39    -0.30       1.00 

30  -0.0k       0.05  -0.11       0.01       0.63       0,15      O.kl       1.00 

kO  0.30  -0.02  O.lk       0.39       0.3k      0.11      0.62       0,77      1.00 

50  -0.29       0.33  -U.k3    -O.kfi    -0.03      0.09       0.03    -0.60    -0.70       1.00 


INTERVAL  * 

c  ■    .OA 

ave 


KESP      1              2  3  k  5            10            20            30           kO            50 

1  1,00 

2  -0.17       1,90 

3  0.k2  -O.kS  1.00 
k  0.57  -0.16  0.67  1.00 

5  0.52      O.08  0.02  0.51       1.00 

10  O.kO      0.20  -O.lk  0.56      0.52      1.00 

20  -0.68       0.2o  -u.3u  -0.32    -0.21    -0.21       1.00 

30  -U.IO  -0,15  -0,13  -0.59    -0.23    -0,37    -0.09       I. 00 

kO  -O.lk      0,0k  0,38  0.k9    -O.lk       O.Ou       O.bU    -0.65       1.00 

50  0.16       0.31  -0.3k  0.26       0.69       0.5S    -O.OC    -0.03    -0.23       1.00 


RESP 

1 

2 

3 

k 

1      1. 

.00 

2       0. 

.31 

l.OU 

3       0. 

.65 

-0.2U 

1.00 

k       0. 

.09 

-0,05 

0.k3 

1.00 

5       0.77       0.31       0.39      0.09       1.00 
lu       0.59       0.19       O.IC    -0.20      0.70      1.00 
/O      0.12       u.6k    -0.51    -0.33       0.33      O.kO 


INTERVAL  5 
p         •   .11 


1.00 


3U    -0.22    -0.k5       0,3U       0,kl    -0,28    -0.38    -0.80       1.00 

ku      0.25      0.65    -0.35    -0.15       0.52       0.k3       0.76    -O.kl       1.00 

5u    -0.3k       0.2k    -O.la       U.3U    -0.k3    -0.37    -0.12       0.k5      0.06       1.00 


INTERVAL  6 

c  •   -.0« 

ave 


RESP      1  2  3  k  5  10  20  30 

1  1. 00 

2  0.37      1.00 

3  0.00      0.06       1.00 
k      O.lk      0.07      0.05       l.OU 

5       0,3b       0.5b       0.01    -0,kk       1.00 
10   -0.55    -0.3b      0.02      0.19    -0.11      1.00 
20    -0.12       0.01      0.31    -O.ek       0.55      0.20       1.00 
30    -0.59    -O.kk    -u.ll      0.73    -0.79       O.kO    -0.57       1.00 
kO    -0.7b    -0.k6    -0.21    -0.2k    -0.17       0.71      0.28       0.20       1.00 
50    -0.52    -0,5k    -0,!'J       0,28    -0.69       O.kO    -0,55       0 ,  5S       0,33       1,00 


INTERVAL   7 


RESP      1              2              3  k              5            10            20           30           kO           50 

1  1,00 

2  U.k3       1,00 

3  -0.36       0.00       1.00 
k    -0.22      U.12       O.kO  1.00 

5      U.39       U.2C    -0.30  -0.37       1.00 

10    -0.57    -0.10       0.k2  0.59    -0.35       1.00 

20       O.kb       0.17    -0.75  -0.26      0.15    -0.15      1.00 

30    -0.«2    -0.26      0.51  0.52    -O.kS      0.58    -0.56       1.00 

kU    -0.78    -0.2fl       0.33  O.OJi    -0.06      0,3i    -0.k7       0.69       1.00 

50    -0.03    -O.kS    -0.32  -0.62    -0.23    -U.52      0.13   -0.12    -0.07      1,00 


TABLES  B.8  through  B.13 

EER  Correlation  Coefficients  for  the  0.3  to  50  Hz  Bandwidth 

The  above  tables  provide  the  correlation  coefficients  betiv'een  the 
ten  wide  band  response  wavefonns  in  Figure  B.l.     These  response  wave- 
forms were  shifted  to  compensate  for  phase  differences  introduced  by 
the  filtering  process.     The  above  intervals  therefore  correspond  to 
those  segmented  to  narrow-band  EER's. 
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CNTl  1  !              3              k 

1  1.00 

;  0.11  1.00 

s  o.as  -0.11    1.00 

b  0.05  -0.22       0.03       1.00 


INTERWU.    1 
D  -    .07 


5       0.50  O.JO      0.55  0.61  1.00 

10    -0.07  -0.70    -0.01  0.20  -0.U«  1.00 

20      0.66  0.23      0.35  -0.1|2  -0.20  O.U  1.00 

50      0.1i5  -O.GO      0.71  0.20  0.38  0.15  -O.IC      1.00 

liO      0.00  -0.2C    -0.09  0.91  0.35  0.1.6  -0.20    -0.01       1.00 

SO    -0,63  -0.UI4    -O.CC  0.19  -O.Cl  0.72  -0.15    -0.1.1      0.1*0       1.00 


30 


kO 


50 


INTEBVW.  2 
c         •  -.02 


1  2  3  k  5  10 

1.00 

0.23       1.00 
-0.k7    -0.17       I. 00 

0.16   -0.79      0.29      1.00 

-0.58       0.59      0.07  -0.88       1.00 
10    -0.36    -0.90       0.25      0.7k    -0.k9       1.00 
20    -0.06    -0.93    -0.03      0.68    -0.58       0.89       1.00 

50    -0.65       0.18       0.75  -0.07       0.k3    -0.07    -0.k3       I. 00 
kO      0.11    -0.82    -0.17       0.5k    -0.57      0.77       0.92    -0.65       1.00 
50      0.11    -0.9k    -0.01       0.83    -0.78       0.88      0.91    -O.kk      0.90       1.00 


CNTL        1              2  3              k              5 
1       1.00 
!      0.75       1.00 

3      O.kt       0.11  1.00 

k    -0.07      0.55  -0.60       1.00 

5       0.1!    -0.37  -0.01    -0.59       1.00 

10    -0.52    -0.93  0.06    -0.7k       0.5k 


I. 00 


HflERVW.  3 
o  •   -.03 


20  -0.76    -0.72  0.06  -0.22  -0.2k  0.63       1.00 

30  -0.b9       0.15  -0.71  0.90  -0.55  -0.k2       0.1k  1.00 

kO  0.2U      0.67  0.10  0.70  -0.77  -0.7k    -0.03  0.53       1.00 

50  0.00      0.65  -0.15  0.86  -0.8k  -0.79    -0.06  0.71       0.95       1.00 


NTL 

1 

1.00 
-0.1k 
-0.92 

0.57 

2 

1.00 

-0.02 

0.66 

3 

1.00 
-0.73 

k 
1.00 

5 

10 

20 

30           kO 
INTtBvm  4 

50 

O.kl 

0.1k 

-0.18 

0.17 

1.00 

10 

0.25 

-0.56 

-0.23 

-0.12 

-0.k6 

1.00 

20 

-O.kS 

-0.22 

0.38 

-0.  35 

-0.83 

0.62 

1.00 

30 

-0.19 

0.89 

0.09 

0.  53 

0.J7 

-0.52 

-0.55 

1.00 

kO 

-0.37 

0.51 

0.05 

0.55 

-0.66 

-0.11 

0.27 

0.21      1.00 

50 

-0.35 

0.85 

0.22 

0.37 

0.20 

-0.85 

-0.39 

0.80      0.5k 

I. 00 

CNTL         1  2  3 

1  1.00 

2  0.0k  1.00 

3  0.0k  -0.27       1.00 


INTERVAL  5 
D  ■   -.07 


k  0.87      0.52    -0.15      l.CO 

5  0.1k    -0.78    -0.27    -0.26  1.00 

10  0.81       0.11       0.5!       0.73  -0.28  1.00 

20  0.67    -0.3S      0.1k       0.k2  0.!k  0.56       I. 00 

30  -0.27    -O.UO    -0.73    -0.39  0.76  -0.75    -0.03       1.00 

kO  -0.98      0.12    -0.11    -0.77  -0.26  -0.79    -0.75      0.21       1.00 

50  -0.66      0.61    -O.kS    -0.27  -0.k2  -0.65    -0.6k      0.15       0.75       1.00 


INTERVAL  6 


1              2  5  k              5            10           20           30 
1. 00 
-0.81      1.00 

0.7k    -0.70  1.00 

-0.18      0.56  -0.23  1.00 

0.07    -0.12  -0.53  0.07       1.00 

10      0.72    -0.82  0.20  -0.k6      0.58       1.00 

20       0.7C    -0.61  0.89  0.11    -0.27      0.20      1.00 

30      0.38    -0.17  0.59  0.30    -0.29    -0.19      0.80      1.00 

kO      0.51    -0.71  0.77  -0.6S    -0.!6      0.28       0.60      0. k6      1.00 

50    -0.89      0.92  ■■0.k8  0.38    -0.k2    -0.9!    -0.53    -0.17    -O.kS       1.00 


INTERVAL    ? 
p  ■   -.07 


CNTl         1              2              3  k  5            10           20           30           kO           50 

1  1.00 

2  -0.11  1.00 

3  -0.20  -0.10      1.00 
k       0.23  0.08      0.81  1.00 

5      0.30  -0.k3      0.75  0.8!       1.00 

10      0.0k  0.k5    -0.85  -0.78    -0.87       1.00 

20      0.02  -0.77    -0.21  -0.k7       0.05       0.06       1.00 

30  -0.81  -0.0k       0.6k  0.15       0.15    -0.!8       0.09       1.00 

kO      0.75  -0.39    -0.5G  -0.18       0.12       0.22      O.kl    -0.80      1.00 

50  -0.9!  0.30      0.05  -0.30    -C.55      0.19    -0.1k      0.70    -0.80       1.00 


TABLES  B.14  through  B.20 

Control  Correlation  Coefficients 

The  above  tables  provide  the  cor- 
relation coefficients  bet\veen  the 
ten  control  waveforms  in  Figure 
3.4  for  eaoh  135  msec  interval. 
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.00 

.05 

.10 

.15 

.20 

.25 

.30 

.35 

.40 

.45 

.00 

1.00 

.05 

0.90 

0.41 

.10 

0.88 

0.53 

0.52 

.15 

0.79 

0.30 

0.28 

0.22 

.20 

0.71 

0.10 

0.35 

0.27 

0.28 

.25 

0.86 

-0.03 

0.17 

0.03 

0.09 

0.11 

.30 

1.07 

0.04 

0.38 

-0.07 

-0.06 

-0.09 

0.00 

.35 

0.58 

0.10 

0.34 

-0.06 

-0.07 

-0.06 

-0.09 

0.00 

.40 

0.63 

0.07 

0.21 

-0.17 

0.10 

-0.12 

-0.04 

-0.08 

0.16 

.45 

0.57 

0.20 

0.06 

0.00 

0.06 

-0.01 

-0.08 

0.01 

-0.13 

0.15 

.50 

0.51 

0.08 

-0.01 

0.01 

-0.21 

0.12 

-0.10 

0.00 

-0.19 

-0.05 

.55 

0.62 

0.13 

0.22 

-0.02 

0.06 

-0.16 

-0.00 

-0.24 

-0.27 

-0.23 

.60 

0.89 

0.20 

0.32 

-0.09 

0.05 

-0.02 

-0.23 

-0.28 

-0.28 

-0.35 

.65 

0.53 

0.27 

0.12 

-0.04 

-0.12 

-0.08 

0.05 

-0.14 

-0.17 

-0.24 

.70 

0.73 

0.25 

0.18 

0.05 

-0.02 

-0.16 

-0.09 

-0.09 

-0.29 

-0.32 

.75 

0.77 

0.30 

0.13 

-0.01 

0.02 

-0.13 

-0.16 

-0.14 

-0.15 

-0.30 

.80 

0.59 

0.26 

0.28 

0.12 

0.03 

-0.02 

-0.05 

-0.21 

-0.14 

-0.16 

.85 

0.57 

0.20 

0.18 

0.18 

0.02 

-0.07 

-0.10 

-0.19 

-0.12 

-0.12 

.90 

0.79 

0.36 

0.14 

-0.01 

0.08 

-0.01 

-0.01 

-0.09 

-0.03 

-0.13 

TABLE  B.21 
Values  for  Cr(t,  t-x)  in  Figure  4.1 

The  following  values  are  measured  relative  to  C^(0,0)  and  the 
values  of  t  and  x  are  given  in  seconds. 
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